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We develop a field theoretical treatment of a model of interacting Brownian particles. 
We pay particular attention to the requirement of the time reversal invariance and the 
fluctuation-dissipation relationship (FDR). The method used is a modified version of the 
auxiliary field method due originally to Andreanov, Biroli and Lefevre [J. Stat. Mech. 
P07008 (2006)]. We recover the correct diffusion law when the interaction is dropped as well 
as the standard mode coupling equation in the one-loop order calculation for interacting 
Brownian particle systems. 

I. INTRODUCTION 

The only existing successful first-principle theory of structural glass transition [1, 2, 3, 4, 5], the 
mode coupling theory (MCT) [6, 7, 8, 9, 10], is beset with absence of controllable approximation 
characterized by smallness parameter. Some years ago one of us attempted to remedy the situation by 
introducing and working out a dynamical fluid model with a Kac-type long range interaction [11] of 
appropriate form among elements of the reference fluid, which is anticipated to exhibit glassy behavior 
[12]. As is well-known, this model has a smallness parameter which is the inverse force range of the 
Kac potential measured in units of inverse microscopic length scale of the reference fluid. However, 
the difficulty with this work is the inadequacy of the expansion scheme which violated the detailed 
balance originating from the time-reversal (TR) invariance of the model equation. 

Recently a great deal of attention [13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28] is 
being paid to go beyond the standard MCT. This stimulated developments of satisfactory perturbative 
calculational methods. Particularly noteworthy is the work of Andreanov, Biroli and Lefevre (ABL) 
[29] where complications associated with the nonlinear TR transformation of the variable set (namely, 
the density field [p] and its conjugate [p] [30]) are avoided by introducing auxiliary variable set [6], [9] 
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which hnearizes the TR transformation. The theory was apphed to interacting Brownian particle 
system [32]. The theory was also applied to continuum nonlinear fluid model[31] where the variable 
set had to be supplemented with the momentum density and velocity fields and their conjugate fields. 
In this case the theory becomes enormously complicated and will not be considered here. However, 
although ingenious is the whole approach, consequences of the theory worked out so far have yielded 
some unsatisfactory features as follows. 

• The equation for the nonergodicity parameter gives non-trivial result even for non-interacting 
Brownian particle systems. 

• The memory integrals entering the equation for the density-density correlation function are 
ill-behaved. 

In our recent communication [33], we have proposed a new set of auxiliary fields still denoted as 
[6], [6] which are defined slightly differently from ABL. However, consequences are drastically different 
so that the two unsatisfactory features mentioned above now disappear. Here we present a detailed 
account of our short communication paper. The paper is organized as follows. In Section II the 
dynamical density field model of interacting Brownian particle system is introduced which is expressed 
as an action integral containing the density field [p] and its conjugate field [p]. This action is shown to 
be invariant under a certain nonlinear TR transformation. This transformation can be converted into a 
linear one by introducing a conjugate pair of auxiliary fields [9] and [6]. The resulting action integral is 
divided into the Gaussian and non-Gaussian parts, each of which is separately TR-invariant. However, 
certain terms coming from the both parts cancel when summed over. Consideration of this fact is 
essential to recover a simple diffusion law for the nonequilibrium averaged density in non-interacting 
case. In Sections HI and IV we develop a renormalized perturbation theory for interacting cases, and 
recover within one-loop order the standard MOT equation for the density-density correlation function 
by invoking the irreducible memory function approach. Section IV summarizes the paper and gives 
discussion. 
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II. THE DYNAMIC DENSITY FUNCTIONAL THEORY 
A. The dynamic equation for the density fluctuations 

We start with the following Langevin equation for the density field p(r, t) of interacting Brownian 
particles 

dtp{r, t) = V • (p(r, V^^) + ry(r, t) (II.l) 

where the Gaussian thermal noise r/(r, t) has zero mean and variance of the form 

< 7y(r, t)i]{r', t') >= 2TV ■ V (^p{r, t)6{r - r')S{t - i')) (11.2) 

where the Boltzmann constant fcs is set to unity, and T is the temperature of the system. Note 
that the noise correlation depends on the density variable, i.e., the noise is multiplicative. This is 
necessary for the Langevin equation (ILl) to satisfy the detailed balance condition so that the system 
is guaranteed to evolve into the equilibrium state governed by the free energy F[p]. In (ILl), F[p] is 
the free energy density functional which takes the following form: 

F[p] =T j dvp{v) (^\n^ -ij + ^ J dr J dr' 6p{r) U{r - r') Sp{r') (II.3) 

where Sp{r,t) = p{r,t) — po is the density fluctuation around the equilibrium density po- In (II.3) 
the first term is the ideal gas part [34] of the free energy, Fid[p], and the second term the interac- 
tion part of the free energy, Fint[p]. Using Ito calculus. Dean [35] has derived the above nonlinear 
Langevin equation for the (microscopic) density of system of interacting Brownian particles with pair 
potential U{r). Earlier, Kawasaki [36] has also obtained the same form of Langevin equation for 
the coarse-grained density with U{r) replaced by — Tc(r), c(r) being the direct correlation function, 
by adiabatically eliminating the momentum field in the fluctuating hydrodynamic equations [31] of 
the simple dense liquids. For this case, the (effective) free energy density functional (II.3) takes the 
Ramakrishnan-Yussouff (RY) form [37] . For further discussions regarding the nature (and controversy) 
of the dynamic equation (II.1)-(II.3), we refer to Ref. [38]. 

B. The dynamic action and the time- reversal invariance 

We consider the corresponding action integral S [p, p] [39] which governs the stochastic dynamics 
of the coarse-grained density variable. The dynamic action can be derived as follows. The average of 
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a dynamic quantity of density variable, A[p\, should be taken over the thermal-noise driven density 
fluctuations satisfying the dynamic equation (II. 1) and (II-2): 

< A[p] >= J Vp A[p] (5 [dtp{v, t)-V- (pV^) - r?(r, t)] )^ 

= jvpjvp A[p] exp (^j dr j dt ip[dtp - V • (pV^)]) (exp J dr J dt ip{r,t)r]ir,t)^ 

= Jt^pJ T^P^ip] exp(<S[p,p]), 

S[p,p] = J drj dt[ip[dtp-V-(^pV^-^)] -Tp{Vpf] (II 

where the auxiliary fleld p is real and the last term involving quadratic p comes from the average over 
multiplicative thermal noise rj. In the first line of (II.4), employing the Ito calculus makes the Jacobian 
of transformation constant. The dynamic action of this form appearing in the above equation with 
the RY free energy functional was first written down in [32]. 

The TR symmetry of the dynamics should be manifested in the dynamic action. In particular, 
under the TR the two fields p and p should transform in such a way that the dynamic action (II.4) 
remains invariant under these transformations. In order to see this TR invariance of the action, one 
can rearrange the dynamic action (II.4) as 

S[p,p] = J drj dt[ip[dtp-V ■ (^pV^-^)] - TipV -{pVip)} 

= I drldt[ip[dtP + iTV-(^pV(^-p+^^-^))]} (II.5) 

This last form of the action suggests that the dynamic action becomes invariant under the field 
transformation 

P(r,-t) = p{r,t) 

p(r,-t) = -^(M)4^ (11.6) 

The invariance of the action under the transformation (II. 6) can be shown as follows. Using the 
transformation of p in (II. 6) one can rewrite (II. 5) as 

= J dr J dtip{r,t)dtp{r,t) + T J dr J dtip{r,t)V ■ (^p{r,t)Vip{r, -t)^ (II.7) 

With the transformation (II. 6) the first integral in (II. 7) is shown to be TR invariant as 

J dr J dt ip{r,t)dtp{r,t) = — J dr j dt ip{r, —t)dtp{r,t) 

= - J drj dt - pir,t) + ^j^^')dtp{r,t) = J dr J dt ip{r,t)dtp{r,t) + ^ J dtdtF[p] 

= J dr J dt ip{r,t)dtp{r,t) (II.8) 
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where the first Une results from the change of integration variable t —>■ —t, and we used the fact that 
the surface term f dt dtF[p]/T vanishes. The second term in (II. 7) is manifestly TR invariant since 
when the time is reversed the spatial integration by parts twice recovers its original form. This proves 
the invariance of the action (II. 4) under the TR transformation (II. 6). Note that the TR invariance 
holds for the generalized form of dynamic action 

S[p,p] = J drjdt [ip[dtp - V • (^(P)V^)] - TipV ■ (P(p)Vip)} (II.9) 

where V(^p{r,t)) is a local function of the density field. This action corresponds to the dynamic 
equation of the form [23] 

dtp{r,t) = V-(p(p(r,t)))V^^)+r?(r,t), 
<r]{r,t)r]{r',t') > = 2TV ■ V (v{p{r,t))6{r - r')5{t - t')^ (11.10) 

There exists another field-transformation leading to the time-reversal invariance of the action. To 

see this, one can again rearrage the dynamic action (II. 4) as 

S[p,p]= J drj dt[ip{r,t){-iT)(^^dtp{r,t) + V ■ {p{r,t)Vp{r,t))^ -ip{v,t)V- (p(r,t)V^)} 

(11.11) 

This form of the action suggests the following transformation 



or equivalently, 



V • (p(r, t)Vp(r, -t)) = ^dtp{v, t) + V- (p(r, i)Vp(r, t)) (11.12) 



p{v,-t) = p{r,t)+iA{p{r,t)) (11.13) 
with the local function A{p{r, t)) defined as 

V • (p(r, t)VA{p{r, t))) = ^5tp(r, t) (11.14) 
One can rewrite (11.11) using the transformation (11.12) as 

S[p, p] = J drjdt [ip{r, t){-T) (V • {p{r, ^)V^/5(r, -t))) - ip{r, t)V ■ {p{v, t)V^) } (11.15) 
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We have previously noted that the first term is manifestly TR invariant. We thus only need to look 
at the last term in (11.15): 

- J drj dtzp(r,t)V- (p(r,t)V^) = - j dr J ^p(r, -i)V • (p(r, t)V^) 

= - J drj dti^V- (p(r,OVp(r,-t)) = - J drj dti^[^dtp{r,t) + V ■ {p{r,t)V p{r,t)) 

= - J J dtip{r,t)V ■ (^pir,t)V^^ (11.16) 

Therefore the dynamic action (II. 4) is TR invariant under the second type of the transformation (11.12) 
as well. The invariance of the action under the second transformation will also hold for the generalized 
form of the dynamic action (II.9). As will be shown in the next subsection, the FDR is readily derived 
from these TR transformations. 

Note that both (II. 6) and (11.12) are intrinsically nonlinear transformations. The equation (II.6) is 
nonlinear because of the noninteracting contribution Fj^, whereas the transformation (11.12) is non- 
linear due to the 'extra' factor of the density field in our dynamic equation (i.e. the multiplicative 
nature of the Langevin equation) . As discussed in the work of ABL [29] , this nonlincarity is the under- 
lying reason why the FDR, obeyed by the action, is not preserved order by order in the renormalized 
perturbation theory developed for the dynamic action (II.4). 

C. Fluctuation-dissipation relation (FDR) 

The FDR is a hallmark of the equilibrium dynamics, which provides a fundamental relationship 
between the correlation of equilibrium fluctuations and the linear response to external perturbation. 
Here it is shown that the FDR is a direct consequence of the TR symmetry. 

The response function R{r, t; r't') is defined as a link between induced density change A < p(r, t) > 
and an external infinitesimal field /ie(r',i') added to F (not to the Langevin equation): 

A < p(r, t) >= J dr' J dt' R{r, t; r't') heir', t') (11.17) 

The contribution of the external field to the free energy, AF = — J dr J dtdp{r,t)he{r,t), will bring 
the corresponding change in the action, AS, which is given by 

AS = J dr J dtip{r,t)V ■ {p{r , t)V heir , t)) (11.18) 

The induced density change is then given by 

A<p(r,0> = (p(r,t)A5)=i J dr' J di' (p(r, i)V' • (p(r', OV'p(r', i')) )/ie(r', *') (11.19) 
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where the integration by parts was performed twice. We see from (11.17) and (11.19) that the response 
function R{r, t; r't') is given by [16] 



i?(r,t;rY) = i{p{r,t)V ■ (p{v' p{v' X, 



(11.20) 



Note that the form of the response function differs from the conventional form of the response function 
i < p{r,t)V'^p{r',t') > which holds for the Langevin dynamics with additive noise: (11.20) reflects the 
multiplicative nature of the original Langevin equation (II. 1) and (II. 2). 

Now we show that the FDR follows directly from the TR transformations. We use the following 
identity 

SS \ I , , , , / . , .,.„, SF 







/5(r,t)- 



(n h 



+ TV ■[p{r\t')V'p{v',t')]+TV ■[p{v\t')V'p{v'X, 



= [p{v,t)[idt,p{v' ,t') - TV ■ {p{v' ,t')V p{v' -t')) +TV ■{p{v',t')Vp{v',t'))]) (11.21) 
where the TR transformation (II. 6) was used. The equations (11.20) and (11.21) give the FDR [40] 



- 7^dtGpp{v -r',t- t') = -Rir - r', t' - t) + i?(r -r',t- t') 



(11.22) 



where Gpp{v — r',t — t') = (^Sp{r,t)Sp{r' ,t')) is the density correlation function. Since the causality 
requires R{r - r',t' - t) = for t > t', (11.22) leads to the standard form of the FDR 



Rir -r',t- t') = -@it - t')^dtGppir -r',t- t') 



(11.23) 



where Q(t) is the Heaviside step function. The FDR (11.22) is more directly obtained from the second 
transformation (11.12) when the second member of (11.12) (denoting the space and time coordinates 
r' and t', respectively) is multiplied by pir,t) and is taken average. 



D. Linearization of the time-reversal transformation 



From now on, we focus on the first TR transformation (II.6). Since the nonlinearity of the TR 
transformation makes the perturbation expansion inconsistent with the FDR [29], this inconsistency 
would be resolved if the transformation is properly linearized. With the form of the free energy given 
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in (II.3), one can explicitly write {1/T)SF/Sp as 

T dp{r, t) po po 

/(<^p(r, i)) ^ - ^ - ( - 5p(r, t)/po)" (11.24) 

n=2 

where f{5p{r,t)) is the contribution of the non-Gaussian (higher than the quadratic) part of -Fidfp]. 
The equation (11.24) leads to an explicit form of the TR transformation (11.6) as 

P(r,-t) = p{r,t) 

p{r,-t) = -p{r,t) +iK *6p{r,t) +if{dp{r,t)), 
k*5p{v,t) = J dr' K{r-r')6p{r',t) (11.25) 

where the kernel K{r) is defined as K{r) = ((5(r)/po ~^U{r)/T^. Note that the transformation (11.25) 
is nonlinear due to the non-Gaussian nature of F^dlp], the ideal-gas part of the free energy. 



1. The Gaussian approximation 

If one entirely neglects f{6p{r,t)) in (11.24), this is tantamount to approximating Fi^lp] to the 
Gaussian form 

F^M^^J dv{5p(y,t)f (11.26) 

In this case, the transformation (11.25) becomes linear and consequently the FDR would be preserved 
by the perturbation theory order by order. However, when (11.26) is substituted to the original 
equation (II. 1), the dynamic equation generates the following two terms 

V • (p(r, i)V^^) = rvV(r, i) + - V • {5pVp) (11.27) 
V op{r,t)J Po 

While the first term is the linear diffusion term, the term which is expected for the noninteracting 
case, the second nonlinear term gives rise to a spurious contribution, incorrectly yielding a nontrivial 
result even in the absence of particle interaction [16, 29, 41]. In order to obtain the correct behavior 
for the noninteracting system, one really should take into account the full logarithmic form of Fj^lp]. 
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2. Introduction of the auxiliary fields 

Fully incorporating Fid[p\ without making truncation, a natural way to make the transformation 
(11.25) linear is to introduce a new field ^(r, t) defined as 

%,t)^/(Mr,t)) = ^^-^ (11.28) 
T dp po 

Note that the definition (11.28) differs from that of ABL in that whereas in the work of ABL, the 
auxiliary field is defined as the functional derivative of the full free energy with respect to density 

9ABL{r, t) = ^SF/Sp{r, t) = K * 5p(r, t) + f{Sp{r, t)), (11.29) 

(11.28) limits the new variable ^(r, t) to the nonlinear part of the transformation. 

With the nonlinear constraint (11.28), using the first member of (11.24) we obtain the ideal-gas 
contribution to the body force as 

V • (pV^) = TV • (pV (^^ + 6^^= TV^p + yV ■ {6pVp) + poTV^9 + TV ■ {dpV9) (11.30) 

We have seen that the first two terms are the contributions from the Gaussian part of Fi^. On the 
other hand, since due to cancellation of the two nonlinear effects the entire ideal-gas contribution to 
the dynamics should be of pure diffusion 



V • (^pVSFia/Sp) = TV ■ (pVln(p/po)) = TV^p, 



the sum of the last three terms in (11.30) should vanish: 

—V • (5pVp) + poTv^e + TV ■ ispve) = o (11.31) 

po 

As shown in the subsection G, this cancellation is used to obtain the correct dynamic behavior for the 
noninteracting case. 

E. New form of the dynamic action 

With the constraint (11.28), the action (II.4) takes the following new form 

5[p,p] = j dvj di{ip[a,p- V- (pV^) - V- (pV^)] -Tp{Vpf} 



S[p, p, e,e] = j dr j dt jip dtp - TV ■ [p^[y +9^^ - TV - {pVU * 6pj 

- Tp{Vp)^+i9{9-f{6p))} (11.32) 
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where the last term comes from the exponentiation of the delta functional 5 d{r, t) — f{p{r, t) 

We now identify the linear transformation under which the new action S[p, p, 6, 6] becomes invari- 
ant. The transformation of p{r, t) and /5(r, t) are already given in (11.25) with f{6p{r, t)) being replaced 
by 9{r,t). The transform of 9{r,t) must be the same as that of p{r,t) since is a local function of p. 
Thus we only need to know how ^(r, t) should transform. Since as shown in (II.7), the terms involving 
p{r, —t) becomes manifestly TR invariant, we only have to consider the term Va{r, t) = ip{r, t)dtp{r, t), 
together with the term Vb{r,t) = i6{r,t)9{r,t). When time is reversed, Va{r,—t) generates a term 
6{r,t)dtp{r,t). This term should be cancelled by the term generated by the Vb{r,—t). This requires 
that 9 transform as ^(r, — t) = 9{r, t) + idtp{r, t). Thus the linear transformation under which the new 
action (11.32) becomes invariant is given by 

p{r,-t) = p{r,t) 

p{r,-t) = -p{r,t) +ik *5p{r,t) + i9{r,t) 

9{r,-t) = e{r,t) 

9{r,-t) = 9{r,t) + idtpir,t) (11.33) 

It is easy to show that the modulus of the associated transformation matrix O is unity (det O = — 1). 

The new action S[p,p,9,9] can be decomposed into the Gaussian part Sg[p, p,9,9] and the non- 
Gaussian part Sng[p, p,9,9] [44]: 

S[p,p,9,9] = Sg[p,p,9,9]+Sng[p,p,9,9] 

Sg[p, p, 9,9] = jdrj dt {ip[dtp - TV^p - poTV^9 - poV^U * Sp\ - Tpo{Vpf + i^^} 
Sng[p, P,9,9] = j dr j dt^ip^-V ■ {SpVU * 5p^ - ^ V • {5pVp) - TV • ((JpV^)] 

-T5p{Vpf ^i9f{5p)] (11.34) 

The actions Sg[p, p, 9, 9] and Sng[p, p, 9, 9] are shown to be separately invariant under the transformation 
(11.33). This is apparent since the transformation (11.33) is linear in Sp, p, 9, and 9. Though with 
the constraint (11.28) the three underlined terms in (11.34) vanish when summed together, one should 
keep each of them in the explicit calculation of renormalized perturbation theory since their presence 
is crucial for the separate invariance of the actions Sg[p, p, 9, 9] and Sng{p, p, 9, 9] which enables one to 
construct the FDR-preserving renormalized perturbation theory from these actions. Nevertheless we 
explicitly show in the one-loop order that the ultimate effect of these three underlined terms is their 
cancellation. We also note that the linearization of the TR transformation inevitably brings back the 
non-polynomial nonlinearity —i9f{5p) in the action. 
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F. New form of the response function 

The presence of the external infinitesimal field he{r,t) leads to shifting f{5p{r,t)) to f{Sp{r,t)) — 
he{r,t)/T in (11.25). This gives rise to a change in action AS[p, p,6,d] = J dr J dti6{r,t)he{r,t)/T. 
Then 

A < p{r,t) >= (^p{r,t)AS[p,p,e,§]) = dt' ^(^5pir,t)e{r',t'))he{r',t') (11.35) 

Therefore with the new form of the action (11.32), the response function takes the form 



R{r,t;r',t') = ^[6pir,t) e{r' ,t')) (11.36) 
Taking correlation of the last member of (11.33) with i5p{r,t)/T we obtain 

l^(^5p{r,t)e{r',-t')) = ^{5p{v-t)e{v'-t')) = ^{Sp{r,t)eir',t')) - ^dt'C,,ir,t;r' ,t') (11.37) 
which is the FDR (11.22). 



G. Dynamics for the noninteracting case: nonperturbative result 

The noninteracting case {U = 0) is an important guide in deaUng with the 9 field since one has to 
recover the linear diffusion law when U = 0. In the absence of the particle interaction, the full action 
(11.32) reduces to 5id[V'] =S[iIj;U = 0] 

Sid[ip] = J dr J dt \^ip\dtp - TV^p - ppTV^O - ■ {ppVp) - TV • {dpVO) 

- TpoiVpf -TSp{Vpf + m-i§f{Sp)^ (11.38) 

where = {Sp, p, 6, 6). We show that the action Sid[i^] yields the dynamic behavior consistent for the 
noninteracting system. We make use of the following identities 

where 1 = (r,t) and 2 = (0,0). The first identity can be written explicitly as 

= (<^P(2)^^) = - rv2)Gp,(l - 2) + 2TpoV\p{l)5p{2)) + 2T{6p{2)V ■ (<5p(l)V/)(l)) > 

(11.40) 

where we used the fact that the sum of the three underlined terms in (11.38) vanishes (see (11.31)). 
Similarly, using the second identity in (11.39), we obtain 

= ('^^(2)^^) = -zporV^</5(l)5p(2)> + z{9{l)dp{2)) - iT{6p{2)V ■ {6p{l)Vp{l))) (11.41) 
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where cancellation of the underlined terms was not used. Since in (11.41), (/5(l)(5p(2)) = (^6{l)Sp{2)) = 
for t > by causality, we obtain 

((5p(2)V- ((5/9(l)V/5(l))> =0 for t>0 (11.42) 

The eqs. (11.40), (11.42), and causahty lead to 

dtGpp{r,t) = TV^Gppir,t), for t > 0. (11.43) 

This result, the diffusion equation for the density correlation function, is the anticipated dynamic 

behavior for the noninteracting system^ [42] <). 

On the other hand, if one employs the ABL approach ((11.29)), 

1 5F 
Tdp{r,t) 

one obtains the following action 

SABLbP] = j drj dt {ip[dtp - TV • (pV^)] - Tp{Vpf + ie(e-k*5p- f{6pfj } (11.44) 

where the subscript ABL of 9 is dropped for simplicity. In the absence of interaction, K{r) = 6{r)/pQ, 
and hence K * 5p{r,t) = 6p{r,t)/po. Therefore 

SABLAdW = J drjdt [ip[dtp - TV • (pV^)] - Tp{Vpf + /(M) } (n.45) 

Applying the above two identities and additional one (^6p{2)5SABL,id[ip]/S0{l)^ = to SABL,id[fp], we 
obtain for t > 



dtGpp{r, t) = TV'Gppir, t) + poTV\f{5p{v, t)Sp{0, 0)) +TV ■ {^6p{r, t)Ve{r, t)Sp{0, 0)) (11.46) 

The ABL action (11.45) thus does not appear to yield the anticipated dynamic behavior for the 
noninteracting particles. This discrepancy is puzzling since no approximation has been made to reach 
(11.46). In fact, however, a careful treatment as detailed in Appendix A leads to the diffusion law 
(11.43). 

III. THE RENORMALIZED PERTURBATION THEORY: THE LOOP EXPANSION 



We are now ready to develop a rcnormalized perturbation theory [39] for the new form of the action 
(11.34), which preserves the FDR order by order. 
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A. Action 

Four field variables '0a(r,t), a = 1,2,3,4 are defined where 

^i{r,t) = 6p{r,t) = p{r,t)-po, Mr,t) = p{r,t), ^3{r,t) = e{r,t), Mr,t) ^ e{r,t) (III.l) 

We arrange fields ■0a(r,t) in column vector and row vector '^^{j) where j stands for a set of 
indices and space-time variable: j = aj,^j with Xj = (rj,tj) . Thus we have 



/ <5p(x,) \ 
Also introduce the notation 



*''0') = ('/'<.j(xy))'^=(5Kxj) p(xj) e(x_j) 9(x_j)) (m.2) 



We formally write the action integral (11.32) as 

1 °^ 1 

'^[V'] = --Go-i(12)^(l)V'(2) + ;^K(12 • • • n)i^(l)m ■ ■ ■ V^(n) 



n=3 



(IIL3) 



(IIL4) 



where T4(12 • • • n) is the fully symmetrized vertices, and we have used the convention that repeated 
arguments j = aj , Xj imply summation over j = aj and integration over Xj . 

B. The one-particle irreducible (IPI) loop expansion 

1. Generating junctionals and Legendre transform 



We introduce a generating functional defined as 

W[J] = In < e'^^^)^^^) >, <•••>= 



(IIL5) 



where J(l) is the local source field. The functional W[J] is the generating functional for the connected 
(cumulant) correlation functions: averages and pair correlations without source field are obtained from 
functional derivatives of W[J] with respect to the source field J in the limit of zero source 



G(12---n) = 



G(l) ^ 0(1) = 
G(12) = 



6J{l)5J{2) ■ ■ ■ 6J{n 
6W[J] 



dJ{l) 



=<V'(i) >, 



dJ{l)6J{2) 



(<V(1)V(2) >-<^(l)<^(2)), etc. (III.6) 
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where response functions are included by using hatted variables, and will not be explicitly mentioned. 

We define a new generating functional known as the vertex functional r[^] via the Legendre trans- 
form 



r[<^] = w[j] - 



We then obtain 



5T[(I^] 5W[J] 5J{2) 5J{2) 



</)(2) - J(l) = - J(l) 



(III.7) 



(III. 



5(t){l) 5J{2) 5(t>{l) 5(t>{l) 
where the first two terms in the rhs cancels since 5W[J]/5J{2) = <p{2). 

The vertex functions r„(12 • • • n) are defined as the derivatives of the vertex functional r[(/)] with 
respect to (/)'s: 



r„(12---n)_ ^^(1)5^(2) ••• 5,/.(n) 
Let us compute the first few vertex functions. Wc already obtained the first one, 



(III.9) 



The second vertex function is given by the inverse of the propagator: 

r2(i2) = - 





SJ{1) 




-1 


r S^W[J] 1 


-1 






ISJii)! 




[6J{1)SJ{2)\ 





G-^(12) 



In order to obtain the higher-order vertex functions, we use the following general equation 
6G-\l2) = -G-\ll')5G{l'2')G-\2'2) = -r2(ll')(^G(l'2')r2(2'2) 



Now the third vertex function Fa (123) is given by 

6T2{12) 6G-^12) 



r3(i23) = 



Using 



we obtain 



or equivalently 



50(3) (5(^(3) 
SG{12) _ 5G{12) 5J{3') 



d(p{3) dJ{3') S(p{3) 



G(123')r2(3'3), 



r3(123) = -r2(ll')r(22')r2(3'3)G(l'2'3' 



(III. 10) 



(III.ll) 



(III. 12) 



(III.13) 



(III. 14) 



G(123) = -G(ll')G(22')G(33')r3(l'2'3') 



(III. 15) 
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Likewise, the 4-point correlation function is given by 

G(1234) = GGTsGT^GG + (2 sym.terms) - GGGGT^ (III. 16) 

where the indices are suppressed for brevity. Therefore one can obtain the higher order correlation 
functions from the vertex functions, which in turn can be obtained from the vertex functional r[^] by 
functional differentiation with respect to 0(1). 



2. IPI Loop expansion 

One can systematically calculate the vertex functional r[0] via the loop expansion [45, 46, 47, 48]. 
We first rewrite (III.5) as 

W[J] = \\n(cj d[^] exp (^S[^] + J(1)V'(1)) ] ) (111.17) 

where c is the normalization factor: it is determined by the requirement iy[J = 0] =0 which leads to 
= J dip exp (^Slipj/X^ . One can perform a formal expansion of W[J] in powers of A (and can set 
A = 1 at the end of calculation) . This is an expansion in the strength of fluctuations and is called the 
loop expansion. We expand the integrand in (III. 17) around the field tpc which extremizes the action 
Slip] + J(1)V'(1) for a given source field J(l): 

Writing ^{1) = tpdl) + \/Ax(1), we expand S[tp] + J(1)V'(1) as 

S[V'] + J(1)V(1) = 5[Vc + VAx] + ^(1)(^c(1) + VAx(1)) 

\n/2 

= (5[V^,] + J(l)V;,(l))+^-p5W(12...n)x(l)x(2)---x(n) (III.19) 

n=2 

where S'c'^^(12 • • • n) = 6^ S[tpc]/Stpc{^)SilJc{2) ■ ■ ■ Stpdn), and the term linear in x vanishes due to the 
relation (III. 18). 

The first two terms in (III. 19) give the leading order results 

W[J] = Wo[J] = S[iPc] + J(l)V'c(l), 
r[<^] = FoM = Wo[J] - J(l)<^o(l) = S[iPc] = S[(l>] (III.20) 
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The functional Wo[J] contains the connected tree (no loop) diagrams only: it is the generating func- 
tional of the connected tree diagrams. The zeroth order result of r[^] is therefore given by the 'average' 

action S [(!>]. 

Considering the expansion (III. 19) up to the term proportional to A, we have the Gaussian integral 
J cixexp[-i|5(2)(12)|x(l)x(2)] = [det |5f (12)1]"^/' = exp ( - ^Trln |5(2)(12)|) (III.21) 
where the identity IndetM = Tr(lnM) is used for the last equality. This means that 

W[J] = Wo[J] + AT^i[J]+0(A2) 

Wi[J] = -l[lVln|5f (12)|-IVln|5(2)(12)|j=o] (III.22) 

where the last term comes from the normalization. 

With the general form of the action (III.4), S'c^^(12) is given by 

5f (12) = -g-i(12; V'c) = -Go '(12) + T4(123)V'c(3) + ^F.(1234)V^e(3)Vc(4) + • • • (III.23) 

Since from (III. 18) ■i/'c = for J = (excluding any symmetry breaking solution of (III. 18) with 
J = 0), we obtain 

Wi[J] = -iTrln(6;-i(ll';Vc)G'o(l'2)) 

= -^TVln (5(12) - F,(ll'3)V'c(3)Go(l'2) - ^F,(ll'34)V'c(3)V'c(4)Go(l'2) + • • •) 
= ^{f,(11'3)Vc(3)Go(1'1) + ^K(ll'34)^e(3)^c(4)Go(l'l) 

- ^K(ll'2)V'c(2)Go(l'2')K(2'3'4)V'e(4)Go(3'l) + ---} (III.24) 

Thus the trace operation generates a set of IPI [49] one-loop diagrams. This expression leads to 

r[0] = Wo [J] + XWi[J] - J(1)0(1) = S[^c] + ^(l)(V'c(l) - </-(!)) + AWi[J] (III.25) 

where we have not yet included the second order contribution to W[J]. We need to express ^c(l) and 
J(l) = -6S[ipc]/Sipc{^) in terms of field. Noting that 

^e(l) = m + 0{X) 

-J(l) = S«(l;^J^gl|l = 5W(l;<^) + S(2)(12;<^)(V',(2)-<^(2))+0(A2), 
S[^P,] = 5[</)] + 5W(l;,^)(Vc(l)- </)(!)) +0(A2), (III.26) 
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we obtain for the first two terms in (III. 25) 

SitPc] + J(l) (Vc(l) - </'(l)) = S[cl>] + 5(1) (1; 0) (^e(l) - 0(1)) + ^^^'^ (12; <^) {M^) - </>(l)) (^^0(2) - 0(2)) 
- (5«(1; 0) + S(2)(12; 0)(V'c(2) - 0(2))) (^e(l) - 0(1)) + 0(A3) 

= 5[0] - ^5(2) (12; cl>){Ml) - 0(1)) (^^0(2) - 0(2)) + 0(A3) (III.27) 

where the two linear contributions cancel. Therefore, up to the one-loop order, the eq. (III.25) leads 
to 

r[0] = 5[0] + XWi[tPc = 0] = -5[0] - ^IVln (g-\(l>) ■ Go) + 0(A2) (III.28) 

where the last term in the last line of (111.27) was not yet included because of its being O(A^). 

We now come to the two-loop calculation. Considering (III. 17) and (111.19), we have the following 
two-loop contributions for W[J] 

W[J] = Wo[J] + \Wi[J] + \^W2[J] + O(A^) 
W2[J] = ^5^(1234) <x(l)x(2)x(3)x(4)>o 

+ iil-^l' S'i\l2?>)Sf\Am) < x(l)x(2)x(3)x(4)x(5)x(6) >o (III.29) 



2! VS!/ 

where VFo[J] and Wi[J] are given respectively by (III.20) and (III. 24). In (III.29), < • • • >o denotes 
the average over the Gaussian distribution 

/d[x](---)exp(-i|5f (12)|x(l)x(2)) 



/d[x]exp(-i|5(')(12)|x(l)x(2) 



< • • • >o= 

Using the Wick's theorem and symmetry of S'c^)(123) and S'c^)(1234), we obtain 

W2[J] = ^5(4) (1234)^(12)^(34) + (123)5(3) (456)^(14)^(25)^(36) 

o 12 

+ J5f)(123)5f)(456)g(12)a(56)g(34) (111.30) 
8 

The diagrammatic expression for VF2[<^] is shown in Fig. 1. While the first two diagrams in Fig. 1 are 
IPI diagrams, the last one in Fig. 1 is IPR (one-particle reducible) diagram. We will see that the 
Legendre transform to r[0] eliminates this IPR diagram, which makes r[0] diagrammatically simpler 
than W[J]. 

Now going back to (III. 25) and adding VF2[J], the second-order contribution of we have for 

the Legendre transform of W[J\ 

r[0] = s[^,] + J(i)(Vc(i) - 0(1)) + Aw^i[J] + \^W2[J\ (III.31) 
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We have seen in (III. 27) that the first two terms in the rhs. of (III.31) has no 0(A)-contribution. Up 
to the second-order in A, they are then given by 

5[^e] + J(l)(V'c(l) - 0(1)) = S[cl>] - ^S(^\12;4>){M^) - <^(l))(Vc(2) - <^(2)) + 0{X^) (III.32) 

Now we compute ■0c(l) — 0(1)- Note that 

= -A5(12;*):5m^+0(A') (III.33) 

where we used (III.18) and (III. 23). Using 
dWi[J](j)] _ 1 6 



6(f)(2) 2 5(/>(2) 



Trln (g-n</')-Go) 



= 4,(rt')5(l'3)«a?l = -i5(3'3)«"<^^'' 



we obtain 



2 ' ^-v y ^^^2) 2^ ' 5(^(2) 

^^?(3'3)^^^^^^ = ^5(3)(233';</>)g(3'3), (III.34) 



V'c(l) - <^(1) = - ^a(12; 0)5(3) (234; 0)^(34; 0) + 0(A2) (III.35) 
Substitution of (III.35) into (III.32) gives 

+ J(l)(Vc(l)-0(l)) 

= 5[0] - ^S^'\12;<P) •( - ^g(13;</.)5(3) (345; 0)^(45; </))) •(- ^^(26; </>)5(3) (678; 0)^(78; </.)) 

= S[^] + ^5(3) (245; 0)^(45; 0)^(26; 0)5(3) (678; 0)^(78; 0) + 0{X^) (III.36) 
8 

where S(^) (12; 0)^(13; 0) = —(5(23) is used. The last term in (III.36) has the same structure as the 
IPR diagram in (III. 30). There is one additional term of the same structure which comes from AVFi[J] 
in (III.31). Since Wi[J] is a functional of i/^dJ) via J(l) = —6S['il^c]/S''Pc{^)i one expands AWi[J] as 

AW^i[J;Vc] = AWi[0] + A^^(Vc(l)-0(l)) + O(A3) 

= AT^i[0] - —5(3) (123; 0)^(23; 0)^(14; 0)^(3) (45g.^)g^5g.^^ ^^^^3) ^jjj ^^^ 
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where the last equality is obtained by use of (III.34) and (III.35). Finally, we consider the last term 
in (111.31): 

\^W2[JM = >?W2[<t>] + o{\^) 

= ^5(^) (1234; 0)g(12; 0)^(34; </.) + ^S^^^ (123; 4>)S^^) (456; ct>)G{U; <k)G{2b; </>)g(36; 0) 

+ — (123; </.)5(3) (456; <^)g(12; </.)g(56; </.)g(34; (^) + ©(A^) (III.38) 

o 

where (III.30) was used. Adding up (III.36)-(III.38), we obtain up to the two- loop order 

m = S[<t>] - ^Tr In (g-\cl>) • Gq) + T^pMI 
TiPi[<j>] = A2{l5(4)(1234;</))e;(12;0)e;(34;</.) + (123; 0)5(3) (456; 0)e;(14; </,)e;(25; 0)e;(36; 0)} 

(III.39) 

Note that the IPR diagram was eliminated in T[(l)]. Therefore, up to the second order in A, Tipi[4>] 
turns out to be the sum of the two-loop IPI diagrams. It was shown that this feature is indeed the 
general structure: rip/[0] is the sum of all (two and higher loop) IPI-diagrams with propagator Q{<p) 
and vertices dictated by the interaction potential. Shown in Fig. 2 are the diagrams for Fi [4>] up to 
the three-loop order [46, 47, 48, 51]. 



C. The two-particle irreducible (2PI) loop expansion 

Although riP7[(/)] consists of IPI diagrams only, some of the diagrams in rip/[0] are 2PR diagrams 
(the diagrams which are disconnected by cutting two lines). Introduction of the bilocal source field 
K{12) can eliminate all 2PR diagrams in the vertex functional obtained via the double Legendre 
transform [46, 47, 50, 51, 52]. Thus the resulting vertex functional F[0, G] has simpler structure than 
F[0]. A specific example is the three-loop result for Ti[(l)]. As shown in Fig. 2, in the three- loop order 
has six IPI diagrams. Among these diagrams, three (the third, fifth, and sixth) diagrams are 
2PI ones (the diagrams which are not disconnected by cutting the two lines) , and the remaining three 
are 2PR ones. The two-loop diagrams in Fig. 2 arc 2PI diagrams as well. These 2PR diagrams were 
shown [46, 47, 50, 51] to be eliminated in the functional r[(/!>, G] obtained via the double Legendre 
transform. 

The generating functional W[J, K] for the connected correlation functions is defined as 

W[J, K] = ln(c f d[iP] exp (5[V'] + ^(1)^^(1) + ^V'(1)^(12)V(2)) ) (111-40) 
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where the normahzation constant is determined by the condition W[J = 0,K = 0] =0. Averages and 
correlations are generated from VF[J, K] as 

SW[J,K] 

0(1) = < V(l) 



5J{1) ' 



^^f"^' ]< 4^{m2)^{3)m > -7 < V'(i)^(2) X mm > (ni.4i) 



5iv:(12)<5iC(34) 4 ^ 4 

The double Legendre transform is then given by 

m G] = W[J, K] - Jimi) - ^K{12) (g(12) + </.(l)0(2)) (III.42) 

where the source fields J and K are should be eliminated in favor of and G. We obtain the derivatives 

SmC] ^ 6W[J,K]6J{2) 6W[J,K]6K{23) ^^(2) . . . . 
I (5J(2) (50(1) 5i^(23) (50(1) (50(1)*^^^ 

1 (5ii:(23) 



2 (50(1) 



■(g(23) + 0(2)0(3)) - K(12)0(2) = -J(l) - K(12)0(2), 



(5r[0,G] 5W[J,K] 5J{3) 6W[J,K]5K{3'i) 5J{3) 1SK{34) 



+ 5K(34) 5G(12) " 5G(U)^^^^ " 2^^021 1^^'^^ + '^^'^'^^"^ J 



(5G(12) (5J(3) (5G(12) (5K(34) (5G(12) (5^(12)"^' ' 2 (5G(12) 

- ^^(12) = -^i^(12) (III.43) 

As before, we first define W[J, K] as 

W[J, K] = \ln(cJ d[ij] exp (sfV-] + ^(1)^^(1) + ^^'(l)^^ (12)V'(2)) ] ) (III.44) 

We then expand the integrand (III.44) around the minimum field ^c(l) of the exponential for a given 
J and K: 

g|i = -(j(l)+i^(12)^,(2)) (III.45) 
Writing ^'(1) = V'c(l) + VAx(1), we obtain 

S[^P] + J(1)V(1) + ^K(12)V'(1)V'(2) = {S[^Pc] + ^(1)^^0(1) + ^K(12)V'c(l)V'c(2)) 

, 00 \n/2 

+ - (5f (12) + (12)) x(l)x(2) + ^ (12 . . . n)x(l)x(2) • • • x{n) (III.46) 

n=3 
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where the Hnear terms in % are eliminated due to (III.45). We see from (III. 45) and (III.46) that the 
formal structure of the expansion is the same as before, provided that the following replacements are 
made: 

S[^,] ^ ^[V'c] + ^X(12)V'c(l)V'c(2), 

Sf\l2) ^ 5p)(12) + K(12) 

g{l2-^c) = -{Si^\l2) + K{12))-^ (III.47) 

There is no change in the higher order terms in (III.46). Therefore the analysis in the previous section 
still holds, and we have 

W[J, K] = Wo[J, K] + Wi[J, K] +W2[J,K] + --- 
Wo[J,K] = S[i^c] + J{l)M^) + lK{12)M^)M2) 

Wi[J,K] = ~Trln(g-\4^c)-Go) 

W2[J,K] = [i^i") (1234)^(12; V'c)^(34;Vc) + ^5^(123)5^ (456)^(14; V'c)^(25; Vc)^(36; V'c) 
+ isf) (123)5(3) (456)g(12;Vc)a(56;Vc)a(34;^e)l (111.48) 

o J 

where the normalization constant is left out. 

The double Legendre transform is now given by 



+ ••• 



r[^,G] = -J{l)4>{l)-^K{l2)(^XG{12) + (l){l)ct>{2)') +Wo[J,K]+Wi[J,K] + W2[J,K] 

= -J{mi) - ^K(12)(aG(12) + 0(l)</,(2)) 

+ SliJc] + J(l)V'c(l) + ^K{12)M^)M^) + Wi[J,K] + W2[J,K] + ■■■ 

= S[A] + (J(l) + K{12)M^)) (Vc(l) - <^(1)) - ^K(12)(V'c(l) - ^(1)) (^c(2) - m) 

- ^Kil2)G{12) +Wi[J,K] + W2[J,K] + --- (III.49) 

Let us look at the first three terms of the rhs in the last line of (III.49) . Using (III.45) and expanding 
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them around (j), we obtain 



1 



^[V'c] + ( + ^(12)^c(2)) (Vc(l) - </-(!)) - 2^(12) (V'c(l) - <^(1)) (V'c(2) - </.(2)) 

= 5[V'c] - ^|||(V'c(i) - '/'(I)) - \k{i2){ui) - <^(i)) (^^0(2) - m) 

= S[4>]+ 5(1) (1; (V-cCl) - ^(1)) + (12; (Vc(l) - <^(1)) (V'c(2) - -^(2)) + • • • 

- (5(1) (1; 0) + 5(2) (12; </>) (V.e(2) - 0(2)) + • • • ) (V'e(l) - 0(1)) 

- ii^(12)(V'c(l)-0(l))(V'c(2)-0(2)) 



-^M - ^ (^^'^ (12; 0) + ^(12)) (Vc(l) - 0(1)) (^-0(2) - 0(2)) 



+ 



(III.50) 



We see that there is no one-loop contribution from these three terms since the terms hnear in (-00 — 
0) cancel in (III.50); they only have the zero-loop, two-loop and higher contributions for r[0, G]. 
Substituting (III.50) into (111.49), we obtain 



r[0,G] = -S[0]-i(5(2)(12;0)+i^(12))(V'c(l)-0(l))(V'c(2)-0(2)) 



- ^if(12)G(12) + Wi[J, K] + W2[J, K] + 



(III.51) 



In order to simplify (III. 51), we first consider the two-point correlations function. It is given by using 
(III.35) and (III.41) 

A 



0(1) 
G(12) 



V'c(l) 



-G{12- (l))S^^^ (234; 0)^(34; 0) + 0{X^ 



,5V'c(l) ^ A (50(1') 6 



5J{2) 5J{2) 2 5J{2) 50(1') 
g(12;V'c) + ^i^(ll';0)G(l'2) + O(A2) 



where F(11';0) = ^ 
by 



g(12'; 0)5(3) (2'34;0)g(34;0) 

(III.52) 

0(12'; 0)5(3) (2'34; 0)^(34; 0) . From (III.52), the inverse of ^(12; Vc) is given 



a-i(12; V'c) = G-'(12) + ^G-i(ll')F(l'2; 0) + 0{\^) 



(III.53) 



Using (III. 53), we obtain the following expression for the term -|K(12)G(12) in (III.51). 

-^i^(12)G(12) = ^G(12)(5f(12) +g-i(12;Vc)) 
= ^G(12)[5(2)(12;0) + 5(3)(123;0)(^e(3) - 0(3)) + G-i(12) + ^G-i(ll')F(l'2; 0)] +0{X') 



= ^G(12)(5(2)(12;0) + G-\12)) + ^G(12)5(3) (123; 0)(Vc(3) - 0(3)) + ^F(22; 0) + 0(A3) 

(III.54) 



A 



A2 
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We now look at the term Wi[J,K], (III.48). Using (III.53), we have 

Wi[J,K] = -^Trln(g-\i;,)-Go)=-^Ti[ln(G-'-Go)+ln{l + ^F) 

= --Tr In (G-^ ■ Goj - —TiF + 0{\^) (III.55) 

where / is the unit matrix. We find that when (III.54) and (III. 55) are added, the last terms in (III.54) 
and (III. 55) cancel against each other. 

Thus we have obtained the following expression for r[^, G] up to the two loop order 

mG] = 5[</>] - ^TVln (g-^. Go) +^Tv(g- 5(2) (</,) + /) 

- ^ (5(2) (12; </,) + if (12)) (^e(l) - 0(1)) (V'c(2) - m) 

+ ^G(12)5(3) (123; - m) + W2[J, K] + 0{\^) (III.56) 

For the last three terms, the same analysis in the previous section that led to (III.39) is applied, and 
the IPR diagram in 1^2 [<^,-^] is eliminated. Consequently, the final expression for r[0, G] up to the 
two-loops with A = 1 is given by 

r[<^,G] = S[<f>] - ^Trln (g'^ ■ Go) + ^'Iv(g • -S^^^ (</>)+ /) +T2Pi[<f>,G] 

r2Pi['P,G] = (1234; 0)G(12)G(34) + -^5(3) (123; (/))5(3) (456; 0)G(14)G(25)G(36) +0(A3) 

o 12 

(III.57) 

where G appears in place of ^(0) since from (III. 52) G = Q{<p) at the lowest order. We can obtain the 
three-loop calculation to explicitly show that the 2PR diagrams contained in r2p/[</>, G] are eliminated. 
It has also been shown [47, 50] that T2Pi[(f>,G] is the sum of the two-loop and higher 2PI diagrams. 
Figure 3 shows the diagrams for T2Pi[(f>, G] up to the three- loop order. 

In the absence of the source fields J = K = 0, we have (/) = 0, S[(f)] = 0, 5(2) (0) = -GqS and 
ST[(I),G]/SG{12) = 0. Also 5(^)(123;0) etc. reduce to the ordinary vertices ^^(123) etc. Therefore we 
have in the absence of sources 

r[0,G] = iTrln(G-Goi) -^Tr(G-Go^-/) +r2P/[0,G] 
The last line is the Schwinger-Dyson equation which can be cast into the form with the self-energy E 

,5r2P/[o,G] 



G-^(12) = Go^-S(12), S(12) = 2- 



(5G(12) 



E(12) = iK(1234)G(34) + iK(134)K(256)G(35)G(46) (III.59) 
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where the last hne, the one-loop expression for the self-energy S, is obtained by use of (III.57). 
Figure 4 shows the self-energy diagrams up to the two-loop order. These diagrams are obtained from 

the diagrams of T2[(p, G] (Fig. 4) by cutting a line (corresponding to taking a derivative of T2Pi[(p, G] 
with respect to G). 

D. Unperturbed propagator 



Recalling we wrote the Gaussian part of the action as 

Sgii^] = -lG^Hl2m)m = -^*^(l) • G^\12) . ^(2), 

we read off the inverse of the unperturbed propagator from Sg[i/j] in (11.34) as follows. 

/ iDi6{12)+ipoVlU{12) o\ 

iDid{12)+ipo\7lU{12) -2poTVld{12) ipoT\7ld{12) 

iTpoVlS{12) -i6{12) 

-i6{12) 



Go '(12) 



V 



(III.60) 



where Di = {d/di + TVf ) and Di = {-d/di + TVf ). 



E. Vertices 



The non-Gaussian part of the action Sng[ip] given in (11.34) can be written as 

Sngm = ^Ki(123)^a(l)V'6(2)Vc(3) + ^Kf^(1234)Va(l)V'6(2)Vc(3)Vd(4) + 

The fully symmetrized vertices V^^^(123) etc. are given by 



K6c(123) 



(III.61) 



y„i(123) = K6c(123) + Kc&(132) + Hc„(231) + Hac(213) + Vcab{312) + Vcba{S21) (III.62) 
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Prom the explicit expression of Sng, we have the following 5 types of cubic vertices 



^ppp(123) 
^;p*(123) 
^p'p(123) 



V>pp(123) 



V^ppe (123) 
%p(123) 



= T/*"*(123) + V' (123) 
H)Vi • [(5(12)ViC/(13) +(5(13)Vii7(12) 



-— V? (5(12)5(13) 



ngl 



Sp{l)dp{2)6p{3) 
dp{l)6p{2)5e{3) 



2TVi • 



(5(12)ViJ(13) 



(-zT)Vi- (5(12)Vi5(13) 



— = -^5(12)6(13) 

6e{l)5pi2)5pi3) py 



(III.63) 



where (5(12) = S{ri — r2)5{ti — ^2), etc. and the indices int and id stand for particle interaction 
and ideal gas, respectively. Note that the vertices V'^*^^(123), V^^^(123), and V^^^(123) are already 
symmetric under exchange of 2 and 3, whereas yppp(123) and V^p0(123) are not. This fact should 
be kept in mind when we compute the one-loop diagrams in order to avoid double-counting of the 
symmetry factor generated from the diagrams. Also U{12) in (III. 63) actually means 

Uil2) = Uiri-r2)diti-t2). 

Thus the above 3-point vertices are nonzero only at the same time ti = i2 = ts- Finally, we observe 
an interesting relationship between the two vertices ¥-^^{123) and V^^^(123): 

PoTVlV^J123) + Vg^{123) = (III.64) 



which will be useful in computing the one-loop diagrams. 

The quartic and higher order vertices come from the term —i9{l)f{dp{l)) in Sng[i>]- For example, 
the nonvanishing quartic vertex is given by 



^^ppp(1234) 



-2^,5(12)(5(13),5(14) 
Po 



(III.65) 



F. Time-reversal symmetry and FDR for G and E 



The time-reversed varible set, (11.33), now denoted by *(i), is given by 



6p{xj[ 



\ 



-/5(xj) + i6{xj) + iK * Sp{xj) 
^(x,) 
e{^j) + idtpi^j) J 



^ 0{dt) ■ *(i) 



(III.66) 
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which also defines 0{dt). The transformation property of the propagator G under TR with the spatial 
coordinates being suppressed is then given by 



G{t-t') = {nm^it')] 

G{t'-t) = (i>m^{t')) = (I^O{dt)^{t)) (0{dt,)^{t')Y) = 0{dt)-G{t-t')-0'^{K' 
= 0{dt)-G{t-t')-0^{-%] 



(III.67) 



By setting i' = in (III. 67), we have 

G{-t) = 0{dt) ■ G{t) ■ O'X-di) 
As for E, using the fact that S transforms like G"^ we obtain 

E(-o = [o^i-dt)] ■ m ■ o-\dt) = o^idt) ■ m ■ oi-di 



(III.68) 



(III.69) 



where wc have used the fact that O, represent time reversal, O^^, {O'^)^^ represent just another 
time reversal which is obtained by changing signs of t. The transformation matrices appearing in 
(III. 68) and (III.69) are explicitly given by 



o{dt) ^ 



( 1 

iK* 




( 1 

i* K 




o\ 

-1 i 
10 

1 y 
o\ 

-1 i 
1 

1 y 



0-10 
Oil 

\o ly 

/ 1 ik* Q idt \ 
0-1 
i 1 

\o ly 



(III. 70) 



where dt is a differential operator acting to the left, and K * {*K) means the spatial convolution of 
K{t) to whatever comes to the right (left). 
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1. FDR for the propagator G 

Straightforward algebra gives the expressions for the elements of the matrices on both sides of 
(III.68): 

Gpp{r,-t) = Gpp{v,t) 

Gpp{r,-t) = ik *Gppir,t) - Gppir,t) +iGp0{r,t) 

Gp0{r,-t) = Gp0{r,t) 

Gp§i^^-i) = -idtGpp{r,t) + G^g{r,t) 

Gpp{r,-t) = ik*Gpp{r,t)-Gpp{r,t)+iGgp{r,t) 

Gpp{r,-t) = ik*(ik*Gpp{v,t)-Gpp{r,t) + iGpe{v,t))-i(^k*Gpp{v,t) + Gpe{v,t^^ 

+ i(ik*Gep{v,t)-Gep{v,t) + iGee{v,t)) =0 

Gpe{r,-t) = ik*Gpe{r,t)-Gpeir,t)+iG0gir,t) 

G-^(r, -t) = ii: * ( - idtGppir, t) + G^^(r, t)) + idtGpp{r, t) + dtGep{r, t) + iG^^iv, i) = 

Gep{r,-t) = Gep{r,t) 

Gepir,-t) = ik*Gep{r,t)-Gep{r,t) + iGe0{r,t) 

Gee{r,-t) = G00{r,t) 

^e§^^^~^) = -i9tGep{r,t) + Ggg{r,t) 

G0pir,-t) = idtGppir,t) + G^^ir,t) 

%(r, -t) = idt {ik * Gpp{r, t) - Gpp{r, t) + iGpe{r, t)) + ik * G^^(r, t) + iG^^{v, t) = Q 

Ggg{v,-t) = idtGpg{r,t) + GgQ{r,t) 

G^-^(r, -t) = idt[- idtGppir, t) + G^^-(r, i)) - i5t%(r, t) = (III.71) 

where we used the property of G^j^{t) = 0. 

Now, various FDRs are obtained from relevant members of (III. 71). First we look at the [p9\ 
element: 

G .{v,-t)-G.{v,t) = -idtGppir, t) 
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Since G^^{r,—t) = for t > by causality we obtain the standard FDR in view of the response 
function (11.36) as 

= ^G^^ir,t) = -@{t)^dtG,,{r,t) (III.72) 

In the same manner we derive the following FDR's for G: 

Gpp{r, t) = iQ{t) (k * Gpp{Y, t) + Gpe{r, t)) 
Gep{r,t) = ie{t)(^k*G0p{r,t) + Gee{r,t)) 

Gg^{r,t) = i@it)dtGep{v,t) (III.73) 

Using the causality requirement Gap{t) = for i > 0, we obtain the FDRs for the adjoint elements: 

Gpp{r,-t) = ieit)(^k*Gppir,t) + Gep{r,t)) =Gpp{r,t) 

= ie{t)dtGpp{r,t) = G^{r,t) 

Gpe{r,-t) = iG{t)(^k*Gpe{r,t) + G0e{r,t)) =G0p{r,t) 

G^^ir, -t) = ie{t)dtGpeir, t) = G^^(r, t) (III.74) 

where Gpe{r,t) = Gop{r,t) was used. Four correlations involving only hatted variables vanish. 
The second line of (III.74) is nothing but the FDR (11.37). 
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2. FDR for S 

As before, we first write out 16 matrix elements of (III.69): 

Spp(r, -t) = ik * Spp(r, t) - idt^poir, t) + iK * (s^^(r, t)+ik* ^pp{r, t) - ^5^S.^(r, t)) 
+ idt (j^gpir, t)+ik* S^-^(r, t) - idt^ggir, t)) = 

Spp(r, -t) = -^pp{r, t) - ik * SppCr, t) - idtT.g.{r, t) 

Epe(r, -t) = iJ^ppir, t) + ik* (iT:^p{y, t) + T.^{v, t)) + idt (^^^--(r, t) + E^^(r, t)) 
^pe(r>-*) = ^pei^,i) + iK*T..^{v,t)+idtT.^^{v,t) 

Spp(r, -t) = -^pp{v, t)-ik* Epp(r, t) + i5tS.^(r, i) 

Ep0{r,-t) = -iT,pp{r,t) -T,p0{r,t) 

Sep(r, -t) = zS^p(r, t) + * (^Spp(r, + Sep(r, t)) - (iS^^(r, t) + S^^(r, t)) 

J:ep{r,-t) = -iT,pp{r,t) -T,gp{r,t) 

^ee{r, -t) = -^pp{r, t) + iT^Opir, t) + iT^poiy, t) = Q 

S,^(r,-t) = zS.^(r,t) + S,^(r,t) 

5^^p(r> -t) = t) + ik* E,-.(r, i) - idtT,^^{v, t) 

^epi^-t) = -^0p{r,t) 
S,-,(r,-t) = zS^-(r,t) + S,-,(r,t) 

S^^(r,-i) = S^^(r,i) (111.75) 
where we used the property Sa/3(r, t) = with unhatted a and /? indices. 
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Using the causality property Sq/jC— = for t > 0, we obtain from (III. 75) the following FDRs: 

^pp{r,t) = ieit)(^ - K *E^p{r,t) + dtE^{r,t)^ 

S^^(r,i) = ze(t)(-i^*S^.(r,t) + atS^,-(r,0) 

S,-,(r,t) = -ie{t)J:^.{r,t) (III.76) 
We obtain the similar FDRs for the adjoint elements using S^^(t) = for i > 0: 
Spp(r,-t) = -ze(t)(x*S^^(r,t)+atS^^(r,t)) 

^ep{r,-t) = -^e(^)E^^(r,^) = Spe(r,t) 

^s§i^,-t) = iem.§{r,t) (III.77) 

Note from (III.25) that while Sp^(r,t) and S^^(r,t) are symmetric under time-reversal, S.^(r,t) 
and E^.(r,t) are anti-symmetric under time-reversal. The unhatted diagonal elements in (III.75) 
vanish. 

G. Dynamical Equation 

The dynamic equations for the correlation and response functions are formally given by the matrix 
Schwinger-Dyson (SD) equation (see (III. 59)): 

Go ^(13) ■ G(32) = (5(12) + S(13) ■ G(32) (III.78) 

where the unperturbed propagator Gq^(12) is given in (III.60). Setting 1 = {r,t), 2 = (0,0), and 
3 = {rs,s), and introducing the space Fourier-transform 

S(r,t) = J S(k,t)e*"^ (III.79) 

where Jj^ = J dk/(27r)^, we can express the above matrix SD equation as 

J ds Go ^(k, t-s)- G(k, s) = 6{t)I + j ds E(k, t - s) ■ G(k, s) (III.80) 
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where I is the 4x4 unit matrix, and the time-integration ranges within (—00,00). The unperturbed 
inverse propagator GQ^(k, t) is calculated from (III. 60) as 

i(dt- poTk'^K{\s.)^S{t) 

- dt - pQTk'^K{k)^ 5{t) 2poTk^S{t) ipoTk'^5{t) 

-ipoTPS{t) 



( 



Go'(k,t) 



V 







Using (III. 81), we write down matrix element of Ihs of (III. 80). 
J dsGQ\k,t-s)-G(k, s) = 



( 



pp 









-i5{t) 



0/ 

(III.81) 



o\ 



X^Gpp + R{2Gpp-iGep) X^Gpp-iRGep X^Gpe + R(2Gpe -iGee) X^G^^-iRG^^ 

-iRGpo - iG^g 

—iGop —iGop —iGee ~'^Gee ) 

(III.82) 



-iRGpp - iG^p 



with 



X+ = i(dt- poTk^K{kfj , X_ = i(^-dt- poTeK{kfj , R = poTk^ 



(III.83) 



We now find the matrix elements of the first term in the rhs of (III. 80). First of all, suppressing the 
wave number and time integration for the moment and denoting a,/3,7 = p, 9, we note the following 

• G]a/3 = ^ajGjf) + YlajGjf) ^ 0, 

whereas 



^■7/3 



since E^^ = and G^^ = by causality. Therefore we get 



-•pp ^pp 



Goo + ^^fl ■ Ga 



pp 



'^pp • Gpp + • G( 



^9p ■ Gpp + T.gg ■ Gq 



■ ^] Op ^ep ■ ^pp + ^ee ' 



^PP ■ Gp0 + E^g ■ Gqq 

'p ■ G]pe 

T.0P ■ Gp0 + Sgg • Ggg 



^PP ■ Gp0 + ^Pd ■ Gg 



^ep ' GpB + ^ee ' ^ee / 
(III.84) 
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where 



L ipp 

[^■Ghp 



^pp ■ Gpp + ^pp ■ Gpp + ^pO ■ Gop + S^g- • 
'^pp • Gp0 + Tjpp ■ Gp0 + • G00 + ■ Gqq 

^ep ■ ^pp + ^ep ■ ^pp + ^96 ■ '^^'p + ^ee ' ^dp 
^§0 ■ ^p^ + ■ + ^ee ■ + ^ee ' ^§e 



(III.85) 



Using the equations of various matrices cxphcitly displayed above we hst below the relevant equations 
of motion where the matrix element is shown as [a/3] 



[21] 



ds 



- dt - poTk^K{k)j Gpp{k, t) + 2poTk^Gppik, t) - ipok^GopiK t) 
^pp{k, t - s)Gpp{k, s) + ^pp{k, t - s)Gpp{k, s) 



+ Epe{k, t - s)Gep{k, s) + S^g(k, t - s)G^^{k, s 



[22] 



[23] 



ds 



ds 



- dt - poTk^Kik)) Gppik, t) - ipoTk^Gepik, t) 
J^ppik, t - s)Gppik, s) + Spe(k, t - s)Ggpik, s)] = Sit) 

- dt - poTk^Kik)) Gpe{k, t) + 2poTk^Gp0{k, t) - ipoTk^GeeiK t) 
^pp{k,t- s)Gpg{k, s) + ^pp{k,t- s)Gpg{k,s) 



+ T.pg{k,t- s)Ggg(k,s) + S^^(k,f - s)Ggg{k,s) 
[24] i(^-dt- poTk^K{k)^ G^{k, t) - ipoTk^G^^ik, t) 

-I 



ds 



^pp{k, t - s)G g{k, s) + Epg{k, t - s)Gg§{k, s) 



= 



iGgp{k, t)- I ds S^_(k, t - s)Gpp{k, s) + S^.(k, t - s)Gpp{k, s 



[41] 

+ S^g(k, t - s)Ggp{k, s) + Eg-^(k, t - s)G§^{k, s) 

[42] -iGgp{k, t)- Jds [S^^(k, t - s)Gpp{k, s) + S^g(k, t - s)Gepik, s) 

[43] -iGgeik, t)- Jds [s^^(k, t - s)Gpe{k, s) + Sg--(k, t - s)Gpe{k, s) 

+ S A„ (k, t - s) Gee (k, s) + (k, t - s) Gnn (k, s) 



[44] 



-iG, 



',{k,t)- I 



ds 



E,(k,t-s)G.{k,s) + ^Uk,t-s)G,g{k,s) =6{t) (111.86) 



We consider [21] element in (III.86) for t > 0. Since Gpp{k,t) = for t > by causality, one can 
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rewrite [21] element as 

dtGpp{k,t) = -poTk''(^K{k)Gpp{k,t) + G0p{k,t)') 

Spp(k, t - s)Gpp(k, s) + Sp0(k, t - s)G0p(k, s) 

Spp(k, t - s)Gpp{k, s) + S-^(k, t - s)G^^{k, s) 



t 

+ i I ds 

' — oo 



+ 



J — oc 



ds 



(III.87) 



where we set the upper limits of the integrals using the causality properties of the self-energies (response 
functions) in the first (last) two integrations: ^^/^(k, t — s) = for s > t and Gap(k.,s) = for s > 
where a,P = p,9. One can simplify the integrals in (III.87) using the FDRs for the self-energies. For 
the first integral, using (111.76) we obtain 

i f ds Spp(k, t - s)Gpp{\i, s) = f ds (if(k)Spp(k, t-s)- 5tS -^(k, t - s)) Gpp{\i, s) (III.88) 

J —oo J —oo 

The second term in the rhs of (III.88) can be rewritten by integration by parts as 

- f ds dt^^ik, t - s)Gpp{k, s) = f ds .^(k, t - s)] Gpp{k, s) 

J —oo J —oo 

= S^g(k, 0)Gpp{k, t) - i;^^(k, oo)Gpp(k, -oo) - / ds ^^^(k, t - s)dsGpp{k, s) 

J —oo 

= - f dsJ:.^{k,t-s)dsGpp{k,s) (111.89) 

J —oo 

where the boundary terms vanish since I]-^(k, 0) = due to the oddness of S-^(k, t) in time (see the 
/5^-element in (111.75)), and Gpp(k, — oo) = 0. The second integral in (III. 87) can be simplified by use 
of (III.76) as 

if dsJ:pe{k,t-s)G0p{k,s)= [ dsE^^(k,t-s)Gep(k,s) (111.90) 

J —oo J —oo 

Therefore we obtain for the first two integrals in (111.87) 

i I ds Y.f,p{k,t- s)Gpp{k,s) + T,pe{k,t- s)G0p{k,s) 

J —oo 

= f ds [Spp(k, t - s) (K{k)Gpp{k, s) + G0p{k, s)) - S-^-(k, t - s)dsGpp{k, s)] (III.91) 

J — oo 

We now show that (— oo, 0)-contributions of (111.91) cancel the last two integrals in (111.87). The third 
integral in (III. 87) involves the response function Gpp{s) with s < 0. Using the p/9-element of (111.71), 
we have 



Gpp{k,-t) = i(K{k)Gpp{k,t) + Ggp{k,t)y for t>0 



(111.92) 
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Setting t = -s{s < 0) in (III.92), we obtain 

Gpp{k, s) = i (JC (k)Gpp(k, -s) + Ggp{k, -s)) , for s<0 

= i(^K{lc)Gpp{k, s) + GepiK s)) , for s < (III.93) 

where we use the fact that the p and 9 field do not change under time reversal. Substituting (III. 93) 
into the third integral in (III. 87), we obtain 

i f ds S^^(k, t - s)G^^(k, s) = - f ds S^^(k, t - s) (if(k)Gpp(k, s) + Gep(k, s)) (III.94) 

J —CO J—QO 

This cancels the (— oo,0) parts of the first two integrals in (III. 91). Similarly, for the last integral in 
(III. 87), we use the FDR, the Fourier transform of the [Op] element of (III. 71) 

G^^(k, -t) = idtGpp{k, t) for i > (III.95) 

Again setting t = —s m. (III. 95), we have 

G^-^(k, s) = -idsGpp{k, -s) = -idsGpp{k, s) (III.96) 

Substituting (III. 96) into the last integral in (III.87), we obtain 

if dsS.^(k,t-s)G^^(k,s)= f cisS.^(k,t-s)5,Gpp(k,s), (III.97) 

J —oo J —oo 

which cancels the (— oo,0) part of the last integral in (III.91). Using these cancellations, we can write 
down the final form for the dynamic equation for Gpp{\s.,t) for t > as 

dtGpp{k, t) = -poTk^ (K(k)Gpp(k, t) + GepiK i)) 
+ 1^ ds [j^ppik, t - s) (K{k)Gpp{k, s) + G0p{k, s)) - S-^(k, t - s)dsGpp{k, s)] (III.98) 

In similar ways, one can obtain dynamic equations for the remaining elements of G in (III. 86). The 
eq. (III. 98) is coupled to the eq. for G0p{k, t) which is given by 

GepiK t) = S^-^(k, 0)Gpp{k, t) + ds S^-(k, t - s) (K(k)Gpp(k, s) + G0p{k, s)) 

- [ dsE^§{k,t-s)dsGpp{k,s) (III.99) 
Jo 

Likewise, we can obtain the two coupled eqs. for Gpo{k,t) and Goe{k,t): 

dtGpe{k,t) = -poTk^(^K{k)Gp0{k,t) + Gee{k,t)^ + ds^pp{k,t - s){K{k)Gpe{k,s) + Gee{k,s)) 

- j dsJ:0{k,t - s)dsGp0{k,s), 

Gee{k, t) = S^^(k, 0)Gpe{k, t) + ds S^--(k, t - s) (K{k)Gpe{k, s) + Geeik, s)) 

- / dsJ:gg{k,t-s)dsGpe{k,s) (Ill.lOO) 
Jo 
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In the same manner, we obtain two sets of coupled dynamic eqs. for the response functions: 

dtGpp{k,t) = -poTk''K{k)Gpp{k,t) - poT^GepiKt) 

+ i 1^ ds [^pp{k,t-s)(^K''{k)Gpp{k,s) + 2K{k)Gpe{Ks) + Ge0{Ks)) 

- 5tS-^(k,t - s)(^K{k)Gpp{k,s) + GpeiKs))' 

GepiK t) = i ds [s,~^(k, t - s) (k'' {k)G pp{k, s) + 2K(k)Gpe(k, s) + Gee(k, s)) 

- 9tS^^(k, t - s) (K(k)Gpp(k, s) + Gpe(k, s)) 
dtG^{k,t) = -poTeK{k)G^{k,t)-poTeGg^{k,t) 

S^p(k,i - s)[K{k)dsGpp{k,s) + dsGep{Ks)) - dtT.^{k,t- s)dsGpp{k,s) 
G^^(k, t) = iS{t) + i 1^ ds [E,^.(k, t - s) (^K{k)dsGpp{k, s) + dsGepiK s)) 



+ i ds 

10 



- dtT,^{k, t - s)dsGpp{k, s) 



(III.lOl) 



H. Static input 



In the present work we are considering the situation in which the system (from some arbitrary initial 
state) has already evolved into the equilibrium state excluding, however, the crystalline state here and 
after. That is, the system is in equilibrium state for t > 0, and the present theory aims to describe the 

dynamics of the equilibrium fluctuations given the static information as input. Therefore the static 
informations of the system enter through the initial values of the dynamic correlation functions: 



Gpp{k, 0) = Sppik), GepiK 0) = Sep{k), G,,(k, 0) = See{k) 



(III. 102) 



where Spp{k) etc. are the equilibrium correlation functions. 

In Appendix A, we derived the following relation between the equilibrium correlation functions: 



K{k)Spp{k) + S0p{k) = l. 



(III.103) 



We also obtain the following relations by setting i = in the dynamic equations (III.99) and (III. 100) 



Ggpik,0) = S^^(k,0)Gpp(k,0) 
G0e{k,O) = S^^(k,O)Gp0(k,O) 



(III. 104) 
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Using the static input (III. 102) and the equihbrium relation (III. 103), one can express S^^(k, 0) and 
G0e(k, 0) from (III. 104) in terms of the static structure factor 5'pp(k) = 5(k) as 

G„{k, 0) = S„(k) = - A'(k))'S(k) (III.105) 

In Appendix A, we showed that SpoCk) = S'ep(k) = in the absence of interaction, U = 0. This 
means that the initial value of the self-energy i;^^(k, t) should vanish for the noninteracting case: 

i;^^(k,0)=0 for U = (III.106) 

If one uses the RY free energy (II. 3), then as shown in Appendix A, K{\s.) is equal to the inverse of 
the static structure factor: i^(k) = S'^^(k). For this case, it implies from (III. 104) and (III. 105) that 
S^^(k, 0) = G^)p(k, 0) = Ge6»(k, 0) = even for the interacting case. 



I. Further nonperturbative results 



We first note that defining 



X{k,t) = K{k)Gpp{k,t) + G0p{k,t) 
Y{k,t) = K{k)Gp0{k,t)+Ge0{Kt), 



(III.107) 



one can rewrite the sets of eqs. (III.98), (III.99) and (III. 100) in terms of Gpp(k,t) and X(k, t), and 
Gp6)(k, t) and Y{k,t) respectively as 

dtGppik, t) = -poTk^X{k, t)+ j ds \^pp{k, t - s)X{k, s) - E^^(k, t - s)dsGpp{k, s)] , 



X(k,i) 



S^-.(k, t - s)X{k, s) - S^^-(k, t - s)dsGpp{k, s) (III.108) 



and 



dtGpoik, t) = -poTevik, t) + ds [Spp(k, t - s)Y{k, s) - S.^(k, t - s)dsGpeik, s) 
Y{k, t) = -g^^Gpeik, t) + £ ds [S^.(k, t - s)Y{k, s) - S^^(k, t - s)dsGpe{k, s)] (III.109) 

In (III. 108) and (III. 109), the first member of (III. 105) was used. One can notice that these two sets 
of equations share the same structure, which leads to interesting results, as shown below. Defining 
the Laplace transform 



/•oo 

G^{k,z)= / dte-'^Gpp{k,t), etc., 

JQ 



(III.llO) 
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we readily obtain the Laplace transform Gpp(k,z) from (III.108) as 

PoTk^ 1 



G^,(k,z) = S(k) 



(l + S^-(k,z))(l-Sf,(k, z)) 

Likewise from another set (III. 109) we obtain 

Gj,(k, .) = S^(k) ■ + . (III.112) 

The two equations (IILlll) and (IIL112) imply that the two correlation functions Gpp(k, and 
Gp6)(k, turn out to be the same when each is normalized with its own initial value (i.e., equilibrium 
correlation): 

GM^t) = . G,,{Kt) = ( Jj-y - i^(k))G,,(k,t) (in.ll3) 

where (IIL103) is used. This relation in turn leads to 



X(k, t) = Kik)Gpp{k, t) + GepiK t) = K{k)Gppik, t) + G^eiK t) = g^G^pih, t) (III.114) 



where the time-reversal symmetric relation ^^^(k, t) = Gp0(k,t) is used. In similar manner, one 
obtains 

y(k,t) = Kik)GpeiKt) + G00ik,t) = ^Gpg{k,t) = ^[-^-K{k))Gpp{^^^^ 
GeeiKt) = - K{k))G,eiKt) = - i^(k))'G,,(k,i) (IILllS) 

The relations (III. 114) and (111.115) with t = are of course fully consistent with (III. 105). Note also 
that (111.113) implies that 

Gp0ik,t) = O for U = (III.116) 

since the relation -fC(k) = 1/S'(k) holds for the noninteracting case. Due to the time reversal symmetric 
relation, we also have 

G0p{k,t) = Gp0(k,t) = Gee{Kt) = O for U = (in.ll7) 



Wc further point out that if the RY free energy is employed from the outset, then wc have from 
(III. 113) and (III. 115) that Gpe(k, t) = Goeik, t) = even for the interacting case since -f^(k) = 1/S'(k). 
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We also note in passing that the second hne of (III. 108) and (III.114) leads to a nonperturbative 
identity 



/' 

Jo 



ds 



5(k) 







(III.118) 



Using (III.114), one can rewrite the dynamic eq. for Gpp{li,t), (III.108) as 



dtCppiKt) = -^^Gpp{k,t) + J\s [j^ppiKt- - ^pei^t- s)dsGpp{k,s)] (III.119) 

Note that the equation now becomes a closed equation for the density correlation function Gpp{k,t) 
alone since the self-encrgics Spp(k, t) and T,^(k,t) can be expressed solely in terms of Gpp(k,t) via 
the FDRs (III. 73) and the relations (III. 113) and the second relation of (III. 115). 



IV. ONE-LOOP REULTS AND MCT EQUATION 
A. One-loop calculations of the self-energies 

We now need to express the self energies appearing in (III. 119) in terms of Gpp(k). Here we 
obtain this expression up to the one-loop order in the loop expansion of the dynamic action. As shown 
in Fig. 4, we have two types of one-loop diagrams. But the first diagram which involves 4-point vertex 
does not contribute since self-energies appearing in (III. 119) are of the form S^^(12), and there is no 
4-point vertex involving two hatted variables. Therefore we only need to consider the second one-loop 
diagrams in Fig. 4, which can be generically written as 

S.^(12) = ^Va^s{lS4)V^^,g,{256)G^y{35)Gs5'{46) (IV.l) 

where indices 7, 7', (5 and 6' include hatted variable indices as well. 

We write down individual expressions of E^^(12) diagrams, starting with S^^(12) 

S^^(12) = lv^J134)V^^^{256)Gpp{35)Gpp{46) (IV.2) 

The diagrammatic expression for (IV.2) is shown in Fig. 5. Note that the symmetry factor generated 
from the exchange 3 <-> 4 and 5 <-> 6 is already contained in the definition of the vertex V^^^. We 
keep this convention for other symmetry-possessing vertices V^^^'*^(123) and V^pp(123). We also note 
that since the vertex V^^^ does not involve U, one-loop expression (IV.2) at t = does not satisfy the 
nonperturbative requirement (III. 106) which should hold for the noninteracting case. Next we write 
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down the expressions for the self-energy S.^(12): 

S«(12) ^ il^;;(134)y^^^(256)G,,(35)G,,(46) 

Eg(12) ^ ^l^,^,(134)y^^^(256)G,p(35)Gpp(46) 

sgV2) = ^yp,e(134)F^^^(256)G,p(35)G,p(46) X 2 

^3^2) ^ ^Fppp(134)y^^^(256)Gpp(35)Gpp(46) x 2 (IV.3) 

We first point out that the diagram E^t^(12) vanishes since Gpp{35) = by causaHty for h = = 
t4 = t and ^2 = *5 = ^6 = with t > 0. Hence we will omit this type of diagrams from now on. The 
multiplication factor 2 in the last two diagrams comes from, e.g., 

Eg)(12) = ^Vppeil34)Vgj256)Gpp{35)Gep{46) + ^Vpep{lS4)V^j256)Gep{35)Gpp{46) 

= \Vppe{lS4)V^j256)Gpp{35)Gep{46) x 2 (IV.4) 

The last line follows from the fact that under the exchanges of dummy variables 3^4 and 5 ^ 6 in 
the second term, the two terms in rhs of the first lines become the same since Vpgp{143) = Vppg{134) 
and V^^^(265) = V^^^(256). The diagrams for non-vanishing S^^(12) is shown in Fig. 6. 
We now move on to the self-energies E^.(12) and Spp(12). 

5^^p(12) = Esg(12) 

i=i 

4^(12) ^ lv^J134)V^fp{256)Gpp{35)Gpp{46) 
Sg(12) = ^V^J134)V^^p{256)Gpp{35)Gpp{46) 
4J(12) ^ ^^^pp(134)F^pe(256)Gpp(35)G,e(46) X 2 

sgV2) = ^y^^^(134)F^,^(256)G,p(35)G,^(46) X 2 (IV.5) 
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The self-energy S^^(12) has 14 diagrams: 



14 



'PP 



12) 



^pp 


:i2) 


= 


-trint 
2 PPP 


[134)V^;,*(256)G,, 


(35)Gpp(46) 




y.(2) 


[12) 


= 


-iT-int 
2 PPP 


[134)Vli^{256)G,, 


(35)Gpp(46) 






[12) 


= 


-trint 
2 PPP 


[134)Vppp{256)Gpp 


(35)Gpp(46) 


X 2 


^PP 


[12) 




2 PPP 


;i34)l^;;(256)G,p 


(35)Gpp(46) 




^PP 


[12) 


= 


2 PPP 


;i34)l^,^,(256)Gp, 


[35)Gpp(46) 




^PP 


[12) 


= 


2 PPP 


[134)Vppp{256)Gpp 


[35)Gp^(46) 


X 2 


^PP 


[12) 


= 


-frint 

2 PPP 


[134)Vpp0{256)Gpp 


(35)Gp0(46) 


X 2 


^PP 


[12) 


= 


2 PPP 


;i34)Fp^e(256)Gpp 


[35)Gp0{46) 


X 2 


^PP 


[12) 




2^pp^ 


;i34)V^;j(256)G,, 


{35)Gep{46) 


X 2 


^PP 


[12) 




2^PP^ 


;i34)l^,V256)G,, 


:35)Gep(46) 


X 2 


^PP 


[12) 




2^pp^ 


;i34)V>pe(256)Gpp 


;35)Ge0(46) 


X 2 


y.(12) 
^PP 


[12) 




\'^ppe 


;i34)V>ep(256)G^e 


;35)Gep(46) 


X 2 


^PP 


[12) 




-j^ppe 


;i34)y^p^(256)Gpp 


:35)Gep(46) 


X 2 


^PP 


[12) 




2^pp0 


;i34)F^^^(256)Gpp 


:35)Gep(46) 


X 2 



The diagrams for Sg^(12) and Spp(12) arc shown respectively in Fig. 7 and Fig. 8. 

It is convenient to classify the entire one-loop diagrams into three classes: the diagrams whose 
vertex biproducts possessing no U (Class 0), the ones possessing single U (Class 1), and the ones 
quadratic in U (Class 2) where we note that U occurs only in V-^^ (the vertex with filled circle in 
Fig. 6 through Fig. 8): 



Class 2 



Class 



Class 1 




(IV.7) 
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First we find it convenient to rewrite the convolution integral in (III. 119) using the identity (III. 118) 



as 



/* 

^0 



ds 



Spp(k,t - ^)^f|^ - S-^(k,t - s)dsG,,{k,s) 



ds [(S^^ - poTk%^) (k, t - ^)%|^ + {poTk^g - S.^) (k, t - s)G,,{k, s)] (IV.8) 
Next, we note the following relation between the two vertices Vg^^{123) and ¥^^^(123) (see (III.64)): 



PoTViV,(123) + Vi^J123)=0 



Opp 



(IV.9) 



This relation reflects the fact that the auxililary field is designed to take care of the nonlinearity of the 
ideal-gas contribution. As shown below, this identity leads to the mutual cancellation of the 5 pairs 
of diagrams in the convolution integral in (III.80). We note the following relation between S^^(12) 
and Til- {12) in the convolution integral: 

/9oTV?S^^(12) + s(J(12) = ^(poTV?y^^^(134) + -l^^^^(134))y^^^(256)Gpp(35)Gpp(46) = (IV.IO) 
The Fourier transform of (IV.IO) can be written as 



s(|(k,i)=PoTfc2s^^(k,i) 



(IV.ll) 



In the same manner wc obtain the same type of relations for the other 4 pairs of one-loop self energies 
as consequences of (IV.9): 



E(;^(k,t) = poTfc^S«(k,t), 45(k,t) = poTk'^f^{k,t), 



45(k,t) = poTk'j^f^iKt), 4?(k,i) = poTfc'4?(k,o, 



pp 

^(8) 

'pp 



(IV. 12) 



^ ' ^' pp^ ' ' Op 

Therefore the kernels Sg^(k, t) and E^^(k, t) are eliminated from the convolution integrals in (IV.8). 
There will be the corresponding cancellations in the higher loop calculations as well. Hence the 
convolution integral (IV.8) can now be rewritten as 

£ ds [Sf,(k, t - ^)%|^ - (S^y + Sg)) (k, t - .)G,,(k, .)] (IV.13) 

where the superscript R denotes the remaining diagrams in the kernel Spp. 

We now compute the remaining diagrams and list the results below. Details of calculation are 
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given in Appendix B. We first write down the results for Class 1 and Class 2 diagrams. 





pI 


h 




-L 


(k 


= 


Po J 


= 


-2iT 1 






'[< 
q 




-^1 


'V 

q 



^2 rr2/ 



r / [(k.q)2c/(q) + (k.q)(k.(k-q))?7(k-q))lG^e(q,i)Gpp(k-q,t) 

(k • q)2 C/(q) + (k • q)(k • (k - q))C/(k - q))! Gep(q, t)Gpp(k - q, t) (IV.14) 



where f = /(i3q/(2^)3. Noting that = S);'^(k,t) (since Gpe(q,t) = ^^^(q,*)), and 



(9), 



PP 



PP 



Gpp(q,t) = i(i^(q)Gpp(q,t) + Gpe(q,i)) with iC(q) = 1/po + U{c^)/T, (FDR (III.73)), we obtain 



2T 



(k • q)2 i7(q) + (k • q)(k • (k - q))C/(k - q))J i^(q)Gpp(q, t)Gpp(k - q, t) 
(k • q)2 t/(q) + (k • q)(k • (k - q))C/(k - q))] Gpp(q, t)Gpp(k - q, t) 
+ 2 / [(k-q)2c/2(q) + (k.q)(k.(k-q))[/(q)C/(k-q))lGpp(q,t)Gpp(k-q,t) (IV.15) 



po .yq 



where the last integral has the same structure as that of S^^^(k, t), (IV.14). The first integral of 
(IV.15) can be simplified by use of 

f (k.q)(k.(k-q))C/(k-q)Gpp(q,i)Gpp(k-q,i)= / (k • q)(k • (k - q))f/(q)Gpp(q, t)Gpp(k - q, i) 

which is obtained by shifting the integration variable q ^ k — q. Using this feature, the first integral 
in (IV.15) is simplified as 

T 
2— 

PO 



(k.q)^ + (k.q)(k.(k-q)) t/(q)Gpp(q, t)Gpp(k-q, i) = 2-fe^ / k-q ?7(q)Gpp(q, t)Gp,(k-q, t) 

PO Jq 

(IV. 16) 

(2), 



which is the same integral as that of EL'^(k, i), (IV.14). Therefore, we obtain 



+ + + E J) + Eg)) (k, i) = ^ k . q C/(q)Gp,(q, t)G,,(k - q, t) 
+ [ [(k-q)2c/2(q) + (k-q)(k-(k-q))i7(q)C/(k-q))lGpp(q,t)Gpp(k-q,t) (IV.17) 

7q 
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Now we calculate and examine Class diagrams. 

Sg(k,t) = / k.qGep(q,t)G^p(k-q,i) 

4?^^'*) = / (k-q)'G,,(q,t)G^^(k-q,t) 

4f = / (k-q)(k-(k-q))Gep(q,OGpe(k-q,0 

^f^'^'^) = / (k-q)'Gep(q,i)Gpp(k-q,i) 

4^'')(k,i) = -2zT2 / (k.q)(k.(k-q))G,p(q,i)G^^(k-q,0 (IV.18) 

7q 

Using the FDR relation, Gep(q, t) = z^K(q)G6)p(q, + Gg0{ci,t)^ , given in (III. 73), we have 

+ = jj^^ • *)^^^(^ - + • q)' i^(q)Gep(q, i)Gpp(k - q, 

= T'Jjk-qf (-^-i^2(q))G,p(q,i)G,,(k-q,i) (IV.19) 
where the last line follows from using (III. 113) and (III. 115): 

2K(q)Gep(q,i) + G,,(q,i) = [2K(q)(-^ - i^(q)) + - K(q))'] G,p(q, i) 



-52(q) 

Likewise using the FDR Gpp{ci,t) = z(K(q)Gpp(q,t) + Gpe(q,t)) , (III. 71), we have 
(4? + 4?) = ^(k • q)(k • (k - q))G,,(q,i)G,e(k - q,t) 
+ /" (k • q)(k • (k - q))i^(k - q)G,p(q, i)Gp^(k - q, t) 

Jq 

= ^(k • q)(k . (k - q)) (-^ - K(q)) • + K{k - q)) Gpp(q, t)G,,(k - q, t) 

- • ^^^"^ • - [(^^ - ^(^)) • (^(k^ + ^('^ - 

+ (^(k^ ~ " '*^) ■ ^ ^^^0] ^^^""^^ " ^' 

= (k • q)(k • (k - q)) _ - ^(q)^(k - q)) Gpp(q, O^PpCk - q, (IV.21) 



where the relation (III. 113) was used. 

Now we are going to add up all the 
care of through (IV. 12). One can recognize the terms sharing with the same structure in the integrand: 



Now we are going to add up all the diagrams of Spp(k, t) where T,^^~ j = 4, 5, 6, 8 are already taken 
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(a) J:f^\k,t) with the first integral of (IV.17), 

(b) {^^/j^^ + with tlic first term in the second integral of (IV.17), 

(c) (S^^^'* + S^^'^^)(k, t) with the second term in the second integral of (IV.17). 
We first compute the following terms: 



(o) 

(i-) 



K(q)A-(k - q; 



T^c(c,)c(k + ^M) + (e(k - q) + 



(IV.22) 



where the direct correlation function c(q) is related to the static structure factor S{q) as S{q) 
/9o/(l — /9oc(q)). Using (IV.22), we obtain 

Epp(k, 0=^2/ [(k • Cifc\q) + (k • q)(k • (k - q))c(q)c(k - q)j G^^(q, ^^^^(k - q, t) 

^2 ^ C/(q)^ 



- 2^ /■(k.q)2(c(q) + ^)G,,(q,t)G,,(k-q,t) 

- ^ /^(k . q)(k . (k - q)) [(c(q) + ^) + (c(k - q) + 
+ / k-q[2C/(q)+Tc(q)lGep(q,t)G^p(k-q,t) 



Gpp(q, t)Gpp(k-q,t) 

(IV.23) 

One can further simplify the last three integrals by noting the third integral becomes by symmetry 



(k.q)(k.(k-q)) c(q) + 



f/(q) 



+ 



(c(k-q) 



+ 



[/(k-q) 



)]Gpp(q,t)Gpp(k-q,i) 



= 2 J (k.q)(k-(k-q))(c(q) + ^)Gpp(q,t)Gpp(k-q,0 



(IV.24) 



Then the sum of the last three integrals is given by 

^the sum of the last three integrals of (IV. 23)^ 



= -2— fc2 f (k.q)(c(q) + ^)Gpp(q,t)Gpp(k-q,0 
+ -k^ [ (k-q)(2;7(q)+rc(q))Gpp(q,i)Gpp(k-q,t) 

PO 7q ^ ^ 

= e / (k-q)c(q)Gpp(q,t)Gpp(k-q,i) 

PO Jq 



(IV.25) 
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One can also rewrite the first integral in (IV.23) as 



(k • q)^c^(q) + (k • q)(k • (k - q))c(q)c(k - q) Gpp(q, t)Gpp(k - q, t) 
^i(k.q)V(q) + i(k.(k-q))V(k-q) 



+ (k-q)(k-(k-q))c(q)c(k-q) Gpp(q, t)Gpp(k - q, 
2 

Using (IV.25) and (IV. 26), we obtain final expression for Epp(k, t) as 



^ ^(k • q)c(q) + (k • (k - q))c(k - q)! ^Gp^iq, t)Gpp{k - q, t) 



(IV.26) 



^ppiKt) = yJ [(k-q)c(q) + (k-(k-q))c(k-q)]'Gpp(q,t)Gpp(k-q,t) 
- —k^ I (k.q)c(q)Gpp(q,i)Gpp(k-q,t) 



(IV.27) 



There are two remaining kernels, T,^}) , (IV.14) and J^'-V , (IV.18) in the convolution integral (IV.13). 
The sum of these kernels is given by 

(E2+S^^))(k,t) = -\ [ k.qC/(q)Gpp(q,t)Gpp(k-q,t)-^ / k • qGep(q, t)Gpp(k - q, t) 
^ P P ^ Po Jq Po Jq. 



.(3) 



^ / k-qc(q)Gpp(q,t)Gpp(k-q,t) 

Po Jci 



(IV.28) 



where (111.113) and C/(q) + ^(^ - K(q)) = -rc(q) is used. 

Summing up, we have obtained the following one-loop expressions of the self-energies in the non- 
perturbative dynamic equation for Gpp{\s.,t) 

dtGpp{k,t) = -P^Gpp{k,t) + J^ds [j^ppiKt - s f'f^^'^ -J:^ik,t - s)Gppik,s) 



S(k) 



^pp(k, t) 



2 
T 



V\k, q) - -V{k, q) Gpp(q, t)Gpp{k - q, t), 
Po 



[ F(k,q)Gpp(q,t)Gpp(k-q,t), 

Jq 



(k • q)c(q) + (k • (k - q))c(k - q) 



(IV.29) 



B. The MCT equation 



The equation (IV.29) takes a different form from the standard MCT equation. While the con- 
volution integral in MCT contains only the time derivative of the density correlation function, the 
first part of the convolution integral of (IV.29) involves the density correlation function itself, instead 
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of its time derivative. Note that this structure is a nonperturbative feature, independent of one-loop 
calculations. It is also very interesting that in the convolution integral of (IV.29) there is an additional 
term which involves the time derivative of the density correlation function. One can appreciate below 
the importance of this term in recovering the standard MOT equation. 

The difference between (IV.29) and the standard MCT is quite analogous to that between the 
reducible and irreducible memory functions [53] appearing in the projection operator approach to 
the dissipative stochastic systems. It has been shown for a general class of dissipative stochastic 
systems obeying the detailed balance condition [54] that the conventional memory function in the 
exact equation for the correlation function, which is obtained from the projection operator method, 
can be further reduced to the so-called irreducible memory function. In particular, the derived exact 
dynamic equation for the correlation function of the slow variable A(t) is given by 

dtCA{t) = -\EA\CA{t)+ f dsMA{t-s)CA{s), 

Jo 

C\{z) = Ca{0) [z + \Ea\ - M^{z)\ (IV.30) 

where jii'yil^^ is a characteristic short time scale in the system, and M^(t) is the conventional memory 
function. The memory function M^(i) turns out to be further reducible to the irreducible memory 
function Mj^(t): 

MA{t) = MY{t)-\EA\-' f dsMA{t-s)MX'{s), 

Jo 

1 + \Ea\-'M^''''''{z) 

(z = [)) can grow 

contrast to M^{z = 0). The above two eqs. lead to the dynamic eq. for C^(i) 



M^{z) = (IV.31) 

Note that M^'*'''"(z = 0) can grow indefinitely when the global relaxation time grows indefinitely in 



dtCA{t) = -\EA\CA{t)-\EA\-^ f dsMY{t-s)CA{s) 

Jo 

C\{z) = Ca{0) [z + l^^l^^^^ ] (IV.32) 

L 1 + \Ea\ ^M/ (zjJ 

For dissipative systems with detailed balance like the one under consideration, the mode coupling 
approximation directly applied to the ususal memory kernel M^(i) in (IV.30) can lead to absurd 
results, which is not the case for the irreducible memory kernel MY{t) in (IV.32). See [54]. This is 
expected to be the same for loop expansions [55]: the dynamic equation for Gpp{\s.,t) with one- loop 
self-energies, (IV.29), is indeed likely to become unstable in the long-time region. 
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Invoking the irreducible memory function formulation, we rewrite (IV. 29) into the correspond- 
ing form of (IV. 32) where the convolution integral involves only the time derivative of the density 
correlation function: 

dtG,p{k, t) = -P^G,,{k, t) - dsM{k, t - s)G,p{k, s) (IV.33) 

The kernel M.{\s.,t), corresponding to the irreducible memory function (A^(k, t) = IEa]"^ MJ'^ (t) with 
\Ea\ = PoTk'^/S{]s.)), satisfies the following equation: 

M{Kt) = t^p{Kt) +^0{Kt) + J^sMiKt - s)ti,p{k,s) (IV.34) 

where I]^^(k, t) = S^^(k, t)/(po7'A;^). It is interesting to recognize the sum of the first two terms in 
(IV.34) is nothing but the standard mode coupling kernel: 



SMc(k,i) = S^^(k,t) + S-,(k,i) 



T 



I [(k-q)c(q) 



+ (k • (k - q))c(k - q) ^Gpp{<i, t)Gpp{\^ - q, t) (IV.35) 



where k = k/A:. 

Now in (IV.34) when M.{k,t) is iterated, the convolution integral generates the terms 
/Q(isSMc(k, t — s)Spp(k, s) + • • • . All these terms are higher-loop terms. We thus see that the 
structure of the theory becomes so simple that only the first contribution A^(k, t) = SMc(k, t) retains 
the one-loop two-particle irreducible structure. Therefore up to the one-loop order, it is perfectly 
legitimate to take this one-loop contribution only, ignoring the terms generated by the convolution 
integral in (III.34). With M{\s.,t) = T^MciKt), the equation (IV.33) reduces to the standard MCT 
equation: 

dtG,,{k,t) = -P^Gp,{k,t) - jys^MciKt- s)Gppik,s) (IV.36) 

V. SUMMARY AND DISCUSSION 
A. Summary 

A field theoretical model for interacting Brownian particle systems is analyzed with a view to 
develop a renormalized perturbation theory consistent with the FDR. This is made particularly trans- 
parent by introduction of a conjugate pair of auxiliary field variables thereby linearizing the TR 
transformation. There is a price to pay for this linearization: a logarithmic nonlinearity reappears in 
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the new dynamical action , the second memeber of (11.32), whereas the original dynamical action, the 
first memeber of (11.32), contains only polynomial nonlinearities. For non-interacting particle systems 
we recover a simple diffusion law as one expects, which is not the case in some recent works [16, 29]. 
For interacting particle cases, we recover the standard MCT in the one-loop order of the renormalized 
perturbation theory. 

B. Discussion 

Having worked out the one loop calculation, the next natural step is to undertake higher loop 
calculations. This is anticipated to be an audacious task: one needs to pay meticulous attention to all 
possible terms (or diagrams) up to desired higher order. Still this is worthwhile since no systematic 
calculation of corrections to the standard MCT is yet available which is crucial to theoretically assess 
successes and inadequacies of the standard MCT, and to push the theory beyond the current MCT. 

In order to systematize higher order calculations, one needs to identify a proper smallness parameter 
(or parameters) of the expansion. Such an expansion parameter is denoted as A in (III. 17), but is 
equated to unity afterwards [56]. 

For a Kac-type system the inverse force range is such a smallness parameter. However, this cannot 
be a whole story as indicated by failure of this approach to meet requirements of TR-invariance 
and FDR that ensues. Non-trivial effects proportional to the absolute temperature arising from 
nonlinearities of the non-interacting part of the action have to be properly included. Development of 
systematic expansion scheme for such model systems would be quite instructive. In future we may 
undertake extensions of our approach to treat multibody correlations [14, 19, 24, 25, 26]. And we 
may also attempt to extend our approach to genuinely non-equilibrium problems like aging [57] and 
colloids under external shear flow [58]. 

Theory of dynamic aspects, not static aspects, is common to old MCT. Its spirit is to predict 
dynamics with statics as input since dynamics is profoundly more difficult. Even this modest aim has 
not been achieved to our satisfaction. 

Now, any satisfactory nonequilibrium theory must describe correctly the dynamics of equilibrium 
fluctuations. Our work is only a first step towards genuine nonequilibrium theory. 

We can compare situation with well- understood critical dynamics. There the dynamical renormal- 
ization group theory (DRG) successfully combines Wilson's RG theory of statics with the old MCT of 
critical dynamics [60] . The resulting DRG is the final theory where indeed statics and dynamics can 
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be treated on equal footing and is capable of predicting even minute details. Extension of the present 
formulation to genuine non-equilibrium situation should be possible. The problems like nonlinear 
rheology are now within our reach. 

It is our hope to develop non-equilibrium theory of glassy dynamics where statics and dynamics 
can be treated on equal footing eventually. However, for that purpose one first needs a theory which 
correctly describes dynamics of equilibrium fluctuations. Sophisticated stage of statics of the current 
liquid state theory [61] gives us a hope to bring dynamics to comparable stage in future. 

Before we end our discussion, it is worthwhile to make comments on the recent related works. 
Szamel reports in a recent paper [28] a diagrammatic formulation of interacting Brownian particles. 
It is a weak coupling expansion scheme (i.e., the expansion in terms of the bare propagator) for the 
hierarchical equations for the multi-density correlation functions. We understand possible usefulness 
of Szamel's approach which does not seem to have problem with the FDR. In some way his (Ander- 
sen's [59] as well) is closer to the original MCT where equilibrium properties are included as input. 
Nonetheless, we wonder how convenient is the calculation. Moreover, some uncontrolled aspects enter 
in obtaining the final hierarchical set of equations. Also the higher order structure is not revealed 
in the formulation. Although the formulation uses 'one-loop approximation', it is not a genuine loop 
expansion theory in contrast to our field theory formulation. 

APPENDIX A: DERIVATION OF THE DIFFUSION EQUATION FROM THE ABL ACTION 
IN THE ABSENCE OF PARTICLE INTERACTION 

We carefully look at the ABL action in the absence of interaction: 

SABL,idm = / / {'P[^tp - ■ {P"^^)] - Tp{yp? + ^^(^ - ^ - f{5p)) } (A.l) 
First note that 

TV • {pve) = V • (pv— ^) = rvV (a.2) 

op 

We rewrite the action (A.l) as 

SablM'^\ = J dr J dt{ip[atP-rvV+rvV-TV- (pVg)]-rp(Vp)^+igfg-^-/(,5p))} (A.3) 

where the two underlined terms cancel due to (A.2). 
We make use of the following identities 

^P^^^ Sp{l) / = °' \^'^^^ 50(1) / = ° ^^-^^ 
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where 1 = {r,t) and 2 = (0,0). The first identity can be written exphcitly as 

= (<5p(2)^^) = z(| - rv2)Gp,(l - 2) + 2rpoV^<p(l)(5p(2)> + 2T{Sp{2)V ■ (<5p(l)Vp(l))> 

(A.5) 

where we used the fact that the sum of the two underUned terms in (A. 3) vanishes. Similarly, using 
the second identity in (A. 4), we obtain 

= ('^/^(2)^^) = -ipoTV^mSp{2)) + z<^(l)5p(2)> - iT{Sp{2)V ■ {5p{l)Vp{l))) (A.6) 

where cancellation of the underlined terms was not used. Since in (A.6), (p(l)5p(2)) = (^(l)(5p(2)) = 
for t > by causality, we obtain 

(5p(2)V- ((5p(l)Vp(l))) =0 for t>0 (A.7) 

The eqs (A.6), (A.7), and causality leads to 

dtGpp{r,t) = TV''Gpp{r,t), for t > 0. (A.8) 

In this way, we recover the diffusion eq. from the ABL action. 
In Sec. II.G, we have obtained ((11.46)) 

dtGpp{r,t) = TV^Gpp{r,t) + poTV^(f{dp{r,t)dp{O,0))+TV • (^5p{r,t)Ve{r,t)5p{0,0)) (A.9) 

Here if the constraint 9 = 5p/po + fi^p) is used, the last two terms cancel: 

poTV^fiSp) + TV • (6pV0) = poTV • (vf + ^VO) = poTV • f - — + — V^) 

^ Po ^ ^ Po Po 

= poTV • (-^Ve - = (A.IO) 
^Po Po > 

where the last equality holds due to (A. 2). Therefore (A.9) reduces to the diffusion eq. (A.8). It 
is extremely puzzling then why ABL's own analysis does not yield (A.8) in the absence of particle 
interaction. 
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APPENDIX B: EQUAL TIME CORRELATIONS 
1. Non- interacting particles 

We begin with the microscopic density p{r) = J2j S{r — Tj). This should not be confused with the 
field conjugate to p. Then 

< p(r)p(rO >= ^ < 6{r - r,)6{r' - r^) >= ^ • • • + J] • • • 

jl 3=1 3^1 

= NS{r - r') < 5{r - ri) > +N{N - 1) < 5{r - ri) >< 5(r - ra) > 

= poS{r-r') + pl (B.l) 

Since < p{r) >= po we get 

Spp{r - r') =< dp{r)dp{r') >= po<5(r - r') (B.2) 
2. Interacting particles 

We use (11.24) 

op(r) Po 

where P = 1/T. We first show that the variable 9 has zero average: 

< > = " < ^ >= " < sS) > " ^ " " s?!^ > 
= -/**^^'"'" = " 

where we have used the fact that dFint/Sp{r) is linear in Sp{r) whose average vanishes. Next we turn 
to another correlation: 

SFid 6p{r') 



S^{r -r')^< Sp{r)e{r') >= {Sp{r) 

= <5(r -r')-K* Spp{r - r') (B.5) 
where K{r) = (^6(r)/po + /3U{r)^ . The Fourier transform of (B.6) is given by 

K{k)Spp{k) + SpeO^) = l 

K{k) = 1 + lt/(k) (B.6) 
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Note that since Spp{ii) = po ((B.2)) in the noninteracting system {U = 0), S'pe(k) vanishes for the 
noninteracting case: 

Spe(k) = 1 - —Po = for U = (B.7) 

Po 

If at the outset the RY free energy functional (II.3) was used, then the function -fC(k) in (B.6) is 
the inverse of the static structure factor Spp(k.) since 

This relation implies from (B.6) that the equal-time correlation Sp0(k) vanishes even in the presence 
of the particle interaction: 

Spe{k) = (B.9) 

Finally we consider 

S„(r - r') = (9(r)e(r')) = ( (/J^ - K . Sp(v)) [l3-^^ - K , 4p(r')) 

= ^"{m^) - 2* •«(■■-■■') + ** ^ • - r') 
= Km^) ^ - r') + . * . S„(r - r') (B.IO) 
This is different from the direct correlation functions as defined by 

jc2 -rex 

where = T — T^^ is the interaction part (excess part) of the full (i.e. renormalized with respect 
to fluctuations) density functional. 



APPENDIX C: CALCULATIONS OF E^^ AND 

Here we derive the equations (IV.14) and (IV. 18). Referring to (IV.3) for E.^ and to (III. 63) for 
the F's we find 
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1. Derivation of (IV. 14) 



1 



^ , (-i)Vi • [5(13)Vit/(14) + 5(14)Vit/(13)] • -S{25)S{26)Gpp{35)Gpp{46) 

^ J 3456 Po 
1 

2Po -^34 



2 , Vi • [<5(13)Vii7(14) + ,5(14)ViC/(13)] Gp^(32)G^p(42) 



= ^ / Vi- [(5(13)Vif/(14)]Gpp(32)G^^(42) 

Po J 34 



Po 

= ^Vi • [Gpp{12)ViU * Gpp(12)] (C.l) 
Po 

where J^^^^q = J d3 J d4 J d5 J d6, etc. With 1 = {r,t) and 2 = (0,0), the spatial Fourier transform 
of (C.l) is given by 

4l^(k,t) = -4 / k-qi7(q)G^^(q,t)Gpp(k-q,0 (C.2) 
which is the first member of (IV. 14). Next we find similarly from (IV.6) and (III.63) 

4^ (12) = l[ ^;;*(134)l^;*(256)Gpp(35)Gp,(46) 

^ J 3456 

= \ I (-OVi • [<5(13)Vii7(14) + ,5(14)Vi[/(13)] 

^ J3456 

• (-i)V2 • [5(25)V2?7(26) + (5(26)V2f/(25)] Gpp(35)Gpp(46) 
= - / Vi[5(13)ViC/(14)]V^(^(25)V^i7(26)] 

J3456 

• [G^^(35)Gpp(46) + G^^(36)G^^(45)] (C.3) 
The first term of (C.3) is computed as 

S(i,-)(12) ^ - / Vi[5(13)ViC/(14)] V^5(25)V^i7(26)] Gpp(35)G,p(46) 

./3456 

= - / V{[Gpp{lb)v{u *Gpp{16)]V\[5{2b)V\U{2Q)\ 

= -V{vi [G,,{12)ViViU * Gpp * U{12)] 

= -ViVi [Gpp{12)V[V{U * Gpp * U{12)] (C.4) 
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where G* implies convolution with the function G(r) as in the case of U*. The second term of (C.3) 
is similarly computed as 



4^^^ (12) ^ - / Vi[<5(13)Vii7(14)] Vi[S{25)ViU{26)] G,,(36)G,,(45) 

J3456 

= - / Vi[Gpp(16)Vi?7*Gpp(15)]V^(5(25)V^?7(26)] 
= - ViV^ [V{U * Gpp{12)V^Gpp * U{12)] 

= -V{v[[v{u*Gpp{12)V[Gpp*Uil2)] (C.5) 

Adding up the above two terms we obtain 

S«(12) = sJ-«)(12) + 4^-^)(12) 
= - ViVi [Gpp{12)V[v{u * Gpp * U{12)] - \7{v[ [ViU * Gpp{12)V[Gpp * C/(12)] (C.6) 

whose Fourier transform is given by 

J:f){k,t) = - [ f(k.q)2c/2(q) + (k.q)(k.(k-q))C/(q)i7(k-q))jGpp(q,t)Gpp(k-q,i) (C.7) 

which is the second member of (IV. 14). Next we have 

Sg(12) = Vif^{134)V^^^{256)Gpp{35)Gpp{46) 

^ 73456 

= 1 f (-i)Vi- [(^(13)Vi?7(14) + (5(14)Vii7(13)] 

^ J3456 

/—iT\ 

X [^r^ [(5(25)5(26)] G,,(35)Gpp(46) 

= Vi-[5(13)ViC/(14)]V^[G,,(32)G,,(42)] 
^ Po^ J34 

= (--)Vi-Vi[Gpp(12)Vi[/*Gpp(12)] 

= ( - -) V?Vi • [Gpp{12)ViU * Gpp(12)] (C.8) 
The Fourier transform of (C.8) is given by 

Eg(k,i) = --k^ [ k.qC/(q)Gpp(q,i)Gpp(k-q,t) (C.9) 
which is the third member of (IV. 14). 
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4?(12) = / ^;,*(134)y^,^(256)G,,(35)G,^(46) 

J3456 

= / (-i)Vi • [5(13)ViC/(14) + 5(14)Vit/(13)] 

J3456 

X (-2r)V2^(25) • [V6<^(56)]Gpp(35)Gpp(46) 

= 2iT [ Vi • [(5(13)ViC/(14)]V2(5(25) • [VeS{56)] 

J3456 

X [Gpp(35)Gpp(46) + Gpp(45)Gpp(36)] (C.IO) 
The first integral of (C.IO) is calculated as 

Sg-")(12) = 2iT [ Vi[<5(13)ViC/(14)] V^J(25)[V[;<5(56)]Gp,(35)Gpp(46) 

^3456 

= -2iT [ Vj[<5(13)V{t/(14)]V^<5(25)Gpp(35)ViGpp(45) 

^345 

= -2iT [ v{[<5(13)V{c/(14)]V^Gpp(32)V^Gp^(42)] 

J 34 

= -2iTV{V'2[Gpp{12)Viv{u * G^^(12)] 

= -2iTViV[ [Gp,{12)V[ViU * Gp^(12)] (C.ll) 
The second integral of (C.IO) is similarly calculated as 

Sj^)(12) = 2iT [ Vi [,5(13)ViC/(14)] V^<5(25)K<5(56)]G,,(45)G,^(36) 

J 3456 

= -2iT [ vi[(5(13)vic/(14)]V^5(25)Gpp(45)V^Gpp(35) 

^345 



= -2iT Vi[,5(13)ViC/(14)]V^Gpp(42)V^Gpp(32)] 

J34 

= -2irV{V^ [ViU * Gpp(12)V^Gpp(12)] 

= -2iTV{v[[v{u * Gpp{12)V[Gpp{12)] (C.12) 

Adding up (C.ll) and (C.12), we obtain 

E(j(12) = s(f)(12) + S2-'')(12) 
= -2iTV{v[ [Gpp(12)ViVif> * Gpp{12)] - 2iTViv{ [V{U * G pp{12)V[G pp{12)] (C.13) 

The Fourier transform of (C.13) is given by 

S|)(k,t) = -2ir / [(k.q)2c/(q) + (k-q)(k-(k-q))C/(k-q))lGpp(q,0Gpp(k-q,0 (C.14) 

This is the 4th member of (IV. 14). 
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4? (12) 



-J ( 

J 3456 



/ 

J3456 



V?;,*(134)F^,e(256)G,,(35)G,,(46) 



(-i)Vi- [(^(13)ViC/(14) + (^(14)Vi;7(13)] • (-ir)V^ (5(25)V^(5(26) Gpp(35)Gp0(46) 
= (-T) / vi[<5(13)V{ [7(14)1 V^5(25)V^(5(26)[Gpp(35)G^e(46) + Gpp(45)G^e(36)l (C.15) 

J3456 L J L L J 

The first integral of (C.15) is computed as 



j:%^\i2) = i-T) 



= i-T) f Vi 



3456 



5{13)V{U{14)\ f(5(25)V^(5(26)l Gpp{35)Gpe (46) 



Gpp{15)V{u * Gp0{16) (5(25)V^(5(26) 

= (-r)ViV^ [Gpp(12)vW{?7 * G^e(12)' 
= (-r)ViVi [Gpp{12)V[v{u * G^e(12)' 

Likewise the second integral of (C.15) is computed as 

5{13)V{U{14)]v'J6{25)ViS{26)\Gpp{45)Gpe{36) 



(C.16) 



4;-'') (12) = (-r) 



3456 



= {-T)V{V[ V[Gpe{12)ViU*Gpp{12) 



Adding up (C.16) and (C.17), we obtain 



(C.17) 



4? (12) 



E?:"Vl2) + E?:*Vl2 



= (-r)ViV'i[G^^(12)V'iVi[/*Gpe(12)J +(-r)VM[V'iGpe(12)ViC/*Gpp(12)J (C.18) 
The Fourier transform of (C.18) is given by 

Sg(k,0 = -r / [(k-q)2c/(q) + (k.q)(k-(k-q))C/(k-q))jG,e(q,t)G^p(k-q,t) (C.19) 
This is the 5th member of (IV. 14). 



Sg)(12) = / F^,,(134)y|;;(256)G,,(35)G,,(46) 

^3456 

T)V{ [<5(13)V'i5(14)j (-i)V2 • [5(25)V2t/(26) + (5(26)ViC/(25)] • Gpp(35)G0p(46) 



3456 



= i-T) [ Vi [,5(13)Vi<5(14)l V'2 k25)V^i7(26)l \Gpp{35)Gep{46) + Gpp(36)Gep(45)l (C.20) 

J3456 L J L J L J 
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The first integral of (C.20) is calculated as 



Sg-")(12) = (-T) 



3456 



(J(13)V{ (5(14)1 [(5(25)V'2i7(26)l Gpp(35)Gep(46) 



= (-T) / Vi Gpp{15)ViGep{16) <5(25)V^i7(26) 

= (-r)ViV'2 [Gpp(12)V^ViG',p * U{12) 

= (-T)ViV'i [Gpp(12)ViViG'9p * C/(12) 

The second integral of (C.20) is calculated as 



J^f^\l2) = (-T) / Vi 5(13)Vi(5(14) (5(25)V^C/(26) G^p(36)Ge^(45) 



(C.21) 



3456 



= (-T) /" Vi[Gpp(16)viGep(15)lv^[(5(25)V^[/(26) 

J56 L -I L 

= (-T)ViV^ [viGep{12)V2Gpp * C/(12) 
= (-r)ViVi \v{G0p{12)V{G,, * C/(12)' 



(C.22) 



Adding up (C.21) and (C.22) leads to 

i:^^}(i2) = e(A-")(12) + S^':'Vl2) 

= (-r)ViV'i [G^^(12)ViViGep * C/(12)l + (-r)vivi [viG0^(12)VlG^^ * [/(12)1 (C.23) 



The Fourier transform of (C.23) is given by 

Eg)(k,t) = -T f [(k.q)2c/(q) + (k.q)(k.(k-q))C/(k-q))]G,p(q,0G^^(k-q,i) (C.24) 
This completes the derivation of (IV. 14). We move on to the derivation of (IV. 18). 



2. Derivation of (IV. 18) 



We start with (12), (IV.4). 



Sg)(12) = / Fppe(134)F^ (256)Gpp(35)G,p(46) 

J3456 

(-ir)vi kl3)Vi(5(14)l • 4<5(25)^(26)Gpp(35)Gep(46) 



3456 '- ^ Po 

^ / Vikl3)vi<5(14)lGpp(32)Gep(42) 

Po -'34 L 
^0 



Gpp(12)ViGep(12) 



(C.25) 
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The Fourier transform of (C.25) is given by 

Sg)(k,t) = -^ J k.qG,p(q,t)Gpp(k-q,0 



(C.26) 



which is the first member of (IV. 18). 

The qauntity T,fP\l2) can be obtained from Ef){12), (C.4), since the vertex Fi^^(256) is obtained 



pp 



from Vi^^{256) with C/(12) being replaced by {T/po)S{12). We thus have 



4f(k>o = f4?(k,t) 



pp 

rp2 



rp2 



PP 

\2 



[/(q)=C/(k-q)=^ 



(k • q)^ + (k • q)(k • (k - q)) Gep(q, t)G,p{k - q, t) 
f k-qGep(q,t)Gpp(k-q,i) 



(C.27) 



which is the second member of (IV. 18). 
Prom (IV.6) 



4? (12) 



3456 



Vppe {lS4)Vppe (256)Gpp {35)Gee (46) 



3456 



(-ir)Vi [<5(13)Vi5(14)l (-ir)V^ k25)V^5(26)l Gpp{35)Geei46) 



= i-T^) / Vi 5(13)V{<5(14) Gpp{32)ViG0e{42) 

J34 L -I L 



= {-T^)V{V[ Gpp{12)V{V[Gge{12 



ivy', 



(C.28) 



whose Fourier transform is given by 



^f}^ (k, t) = -T^ [ (k • af Geeiq, t)Gpp{k - q, t) 



(C.29) 



This is the third member of (IV. 18). 
Again from (IV.6) 



E^f (12) = / Vppg{134)Vpep{256)Gpg{35)Gep{46) 

J 3456 

= [ Vpp0{134)Vpp0{256)Gp0{36)Gep{45) 
J3456 

= (-T2) / v{\s{13)v{S{14)]vl,\5{25)Vi,5{26)\Gpg{36)G0p{45) 
J3456 ^ 



34 



(5(13)vi5(14)l fGep(42)V^Gpe(32) 



{-T^)V{V[ \v{Gep{12)v[Gpe{12) 



(C.30) 
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whose Fourier transform is given by 

4f (k,t) = -T2 [ (k.q)(k.(k-q))G,p(q,t)Gpe(k-q,t) 
which is the 4th member of (IV. 18). 



(C.31) 



S(f (12) 



3456 



3456 



Vpp0{134)Vppp{256)Gpp{35)Gep{46) 

i-iT)v{ kl3)Vi (5(14)1 • (-2r)V^(5(25)V^(5(56)Gp^(35)GeM46) 



'34 



= i-2iT')ViV[ [Gpp{12)ViV[G0p{12) 

The Fourier transform of (C.32) is given by 

E^f (k,i) = -2iT2 / {k.qfG0^iq,t)Gpp{k-q,t) 

which is the 5th member of (IV. 18). 
Finally, 

S(^)(12) = / VppeilS4)Vppp{256)G,^{35)Gep{46) 

J3456 



(C.32) 



(C.33) 



3456 



3456 



{-2iT'^ 



Vpp0{134)Vppp{256)Gpp{36)G0p{45) 



iT)v{\s{13)V{S(14)\ ■ {-2T)\V25{25)Vq6{56)\Gpp{36)G0p{45 



34 



Sil3)v{S{14)\ \G0pi42)ViGpp{32) 



= {-2iT^)ViV[ [v[Gpp{12)V{G0p{12) 
The Fourier transform of (C.34) is given by 

^f^\k,t) = -2iT^ [ (k-q)(k.(k-q))Gep(q,OG'pp(k-q,0 
which is the last member of (IV. 18). This completes the derivation of (IV. 18). 



(C.34) 



(C.35) 
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FIGURE CAPTIONS 

Figure 1. Diagrammatic expression for W^2['^] ((in.30)). The line denotes the propagator Q{tpc)-, and the 
crossing points with 3-branching lines and with 4-branching lines denote respectively the vertices 
si^^ and si'^\ While the first two diagrams are IPX diagrams, the last one is IPR diagram. 

Figure 2. Two and three loop diagrams for rip/[^]. The line denotes the propagator Q{(j)), and the crossing 
points with 3-branching lines and with 4-branching lines denote respectively the vertices S^^^ {(j)) 
and S^^\(f)). The two loop diagrams are both IPI and 2PI diagrams. Other 2PI diagrams are 
the 3rd, 5th, and 6th diagrams. 

Figure 3. Two and three loop diagrams for T2Pi[(t>,G]. Via the double Lcgcndre transform the 2PR 
diagrams in Fig. 2 are eliminated in T2P1 [(p, G] . Here the line denotes the full propagator G. 
In all figures, the vertices with more than 4-legs are not shown. These higher vertices do not 
contribute to the two- loop results for both riP7[^] and r2P7[<^], and hence do not contribute to 
the one-loop result for the self-energy. 

Figure 4. One and two loop diagrams for the self-energy S(12). These diagrams are obtained from those 
in Fig. 3 by cutting a single line. 

Figure 5. The one-loop diagram for S^^(12). 

Figure 6. The one-loop diagrams for S^(12). The filled circle in the vertex ppp denotes the vertex V^^p, 
which is to be distinguished from the vertex V?^^ (without filled circle) 

Figure 7. The one-loop diagrams for E^.(12). 

Figure 8. The one-loop diagrams for S^^(12). 
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